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CHAPTER  I 


INTRODUCTION 


This  dissertation  deals  with  the  use  of  quadri 1 ateral  surface- 
patch  currents  in  the  moment  method  solution  of  antenna  and/or  scatter¬ 
ing  problems.  Wires  and  wire/plate  attachments  are  also  considered 
here.  One  of  the  goals  in  this  work  is  to  develop  a  general  purpose 
computer  code  which  can  solve  near-  as  well  as  far-zone  parameters 
of  arbitrarily  shaped  antennas/scatterers  accurately. 

The  basic  formulation  of  the  problem  is  in  terms  of  an  integral 
equation  where  the  unknowns  are  the  surface  currents.  Here,  Rumsey's 
Reaction  Integral  Equation  (RIE)  Cl3,  which  is  a  general  form  of 
the  electric  field  (EFIE)  and  magnetic  field  (MFIE)  integral  equation 
C23,  is  used.  The  choice  of  electric  test  source  (used  in  this  work) 
reduces  the  RIE  to  the  EFIE  which  is  applicable  to  closed  bodies  as 
well  as  open  bodies,  such  as  infinitely  thin  plates  and  wires.  On 
the  other  hand,  if  magnetic  test  sources  are  used,  then  the  RIE  reduces 
to  the  MFIE.  Although  the  MFIE  is  only  applicable  to  closed  bodies, 
it  is  much  simpler  to  implement. 

Considerable  work  on  the  thin-wire  formulation  resulted  in  many 
useful  and  versatile  computer  codes  [1 3  ,  4,  and  53.  These  codes  are 
ideal  for  thin-wire  problems  and  can  be  used  to  model  solid  surfaces 
by  representing  them  in  terms  of  a  wire-grid.  Richmond  C63  was  the 
first  to  show  the  usefulness  of  the  wire-grid  models  for  scattering 
from  planar  as  well  as  three-dimensional  objects.  Numerous  work  on 
this  subject  followed  as  summarized  inC73.  The  advantage  of  wire- 
grid  modelling  is  its  ability  to  treat  solid  objects  with  sharp  edges, 
corners,  and  thin-wire  appendages.  However,  the  wire-grid  model  has 
to  be  constructed  with  care  in  order  to  have  the  minimum  number  of 
unknown  wire  currents,  to  accurately  represent  the  physical  object; 
yet  not  violate  the  assumptions  of  the  thin-wire  approximation.  Anoth¬ 
er  disadvantage  of  the  wire-grid  modelling  is  that  it  may  not  reliably 
predict  accurate  near-zone  parameters  such  as  surface  currents  and 
impedances. 

Early  work  on  surface-patch  modelling  was  done  using  the  MFIE  to 
solve  closed  surface  scatterers.  0shiroC83  used  a  pulse  basis  func¬ 
tion  and  point  matching  approach  to  solve  the  MFIE  for  various  well  de¬ 
fined  three-dimensional  objects.  This  method  was  extended  for  arbitrar¬ 
ily  shaped  surfaces  by  Knepp  and  Goldhirsh  C 93  through  the  use  of  non- 
planar  cells.  In  their  work,  nine  points  on  the  body  were  used  to 
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find  the  directions  of  the  two  tangential  surface  current  expansions 
for  each  cell.  Albertsen,  etal . ,  [103  allowed  wire  attachments  to 
a  closed  body  by  applying  the  MFlE  for  the  body  and  the  EFIE  with 
the  thin-wire  approximation  for  the  wires.  Planar  quadrilateral  patch¬ 
es  were  used  to  simulate  the  body  along  with  pulse  basis  functions 
with  point  matching  to  represent  the  surface  currents.  Special  consid¬ 
erations  were  also  given  to  the  surface-patches  on  which  wires  are 
attached.  Based  on  this  wire/closed  body  formulation,  the  user-oriented 
Numerical  Electromagnetic  Code  (NEC)  was  developed  by  Burke  and  Poggio 
C113  for  solving  electromagnetic  problems  involved  with  wires/closed 
bodies.  In  this  code,  the  user  specifies  the  shape  of  the  body  by 
first  dividing  it  into  surf ace-patches .  The  area,  x-y-z  coordinates 
of  the  center,  and  the  unit  normal  vector  of  each  patch  are  then  enter¬ 
ed  as  inputs  into  the  computer  code.  Recently,  J.J.H.  Wang  [[123]  used 
the  planar  triangular  patches  along  with  the  MF1E  to  solve  the  scatter¬ 
ing  problem  of  closed  bodies.  The  basis  function  is  related  to  the 
incident  field  by  the  physical-optics  current  technique.  Consequently, 
the  moment  method  impedance  matrix  is  dependent  on  the  incident  field. 

N.N.  Wang,  et  al .,  C 133]  employed  the  RIE  (using  electric  test 
source)  along  with  rectangular  sinusoidal  surf ace-patch  basis  function 
and  Galerkin's  method  to  solve  for  the  currents  and  scattering  paramet¬ 
ers  of  planar  and  non-planar  surfaces.  Sankar  and  Tong  [[143] used 
the  finite  element  approach,  where  the  equation  for  the  surface  currents 
is  made  stationary  with  respect  to  the  set  of  trial  functions  (similar 
to  Galerkin's  method  for  EFIE).  The  planar  triangular  patch  has  been 
used  because  of  formulation  ease  and  better  representation  of  arbitrary 
surfaces.  Planar  triangular  patches  were  also  used  in  the  EFIE  formula¬ 
tion  by  Wilton,  et  al . ,  C153.  Again,  pulse  basis  function  and  point 
matching  were  used  in  the  moment  method  solution.  Newman  and  Pozar 
C16U  modelled  composite  wire  and  surface  geometries  by  using  sinusoidal 
rectangular  surface-patch,  thin-wires,  and  wire/plate  expansion  modes. 
Galerkin's  method  was  used  to  solve  the  RIE  (with  electric  test  source) 
for  far-  as  well  as  near-zone  parameters  such  as  surface  currents 
and  impedances.  The  rectangular  sinusoidal  surface-patch  was  extended 
to  quadrilateral  and  triangular  surface-patches  by  Singh  and  Adams 
C173. 


The  basic  formulation  in  this  work  parallels  much  of  the  work 
done  by  Newman  and  Pozar  H 163,  except  that  quadri 1 ateral  surface-patch¬ 
es  are  now  used.  Attention  is  also  given  to  the  problem  of  automatical¬ 
ly  generating  the  surf ace-patch  modes  on  a  given  polygonal  plate. 

One  of  the  goals  in  this  work  is  to  produce  a  user-oriented  computer 
code  for  electromagnetic  problems  dealing  with  plates  and/or  wires. 

The  general  idea  is  to  let  the  user  model  a  solid  object  by  a  series 
of  connecting  plates,  and  only  the  coordinates  of  the  plate  corners 
are  needed  to  specify  the  plates.  The  automatic  generation  of  surface- 
patch  modes  by  the  code  will  greatly  reduce  the  amount  of  input  data 
required.  It  should  be  noted  that  the  basic  surface-patch  dipole 
can  be  either  a  planar  or  a  non-planar  V-dipole  which  can  be  used 
to  model  any  arbitrary  three-dimensional  surface. 
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Chapter  Two  presents  the  reaction  for/nul ation  for  the  moment 
method  solution  of  composite  surface/wire  geometries.  The  wire,  plates, 
wire/plate  attachment,  and  plate  overlap  modes  are  defined.  Techniques 
for  accurate  and  efficient  evaluations  of  the  mutual  impedances  are 
discussed. 

Chapter  Three  describes  Mie  method  for  automatical ly  generating 
surface-patch  modes  for  a  given  polygonal  plate.  >rit.eria  that 
must  be  satisfied  in  order  to  have  a  set  of  surf ace-patch  modes  which 
satisfy  the  current  conditions  on  the  plates  jre  listed.  In  addition, 
the  geometric  features  of  the  surf ace-patch  which  will  allow  efficient 
evaluation  of  the  impedance  matrix  <e-c  d>>sc  r  dire.  A  standard  subdivi¬ 
sion  technique,  as  well  as  two  other  oigoi  i thins  i or  automatic  surface- 
patch  mode  generation,  are  presented. 

Chapter  Four  shows  the  results  obtained  by  using  the  subdivision 
methods  presented  in  Chapter  Three  for  different  polygonal  plates. 

Both  input  impedance  and  backscatter  data  are  compared  with  other 
independent  calculations  and  measurements. 

Chapter  Five  describes  the  use  of  the  moment  method  to  analyze 
the  microstrip  antenna.  The  microstrip  patch  is  modelled  by  the  sur¬ 
face-patch  mode  and  the  dielectric  is  modelled  by  volume  polarization 
currents.  The  volume  polarization  currents  are  related  to  the  surface 
current  by  the  boundary  condition  at  the  conductor  and  dielectric 
interface.  This  results  in  having  the  surface  currents  as  unknowns 
and  the  dielectric  is  taken  into  account  by  a  modification  of  the 
free-space  impedance  matrix.  Calculations  are  made  for  many  patch 
geometries  and  the  results  are  compared  with  measurements  and  calcula¬ 
tions  by  other  investigators. 

Chapter  Six  is  the  conclusion  and  includes  suggestions  of  related 
topics  for  further  research. 
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CHAPTER  II 


FORMULATION  OF  THE  PROBLEM 


In  this  chapter,  the  electromagnetic  scattering  or  antenna  problem 
is  formulated  in  terms  of  a  surface  integral  equation  based  on  the  reac¬ 
tion  concept  C13  .  The  unknown  is  the  surface  current  which  can  be  cal¬ 
culated  using  the  moment  method.  The  surface  structure  of  the  antenna 
or  scatterer  is  modelled  by  a  combination  of  plates  and  thin  wires. 

Thus,  the  surface  currents  can  be  represented  by  a  finite  series  with 
three  basis  functions:  surface  dipole  modes,  wire  dipole  modes,  and 
wire-to-surface  attachment  modes. 


A.  The  Reaction  Formulation 


The  use  of  the  reaction  concept  to  derive  an  integral  equation 
which  describes  the  electromagnetic  fields  in  an  antenna  or  scattering 
problem  has  been  presented  by  Richmond  093,  Wang,  et  al.,  C 133  and  New- 
mand  and  Pozar  £16],  but  is  included  here  for  completeness.  The  dis¬ 
tinction  between  an  antenna  and  scattering  problem  is  that  when  the 
source  is  located  on  or  near  the  object,  it  is  considered  as  an  antenna 
problem  and  when  the  source  is  distant  from  the  object,  it  is  viewed  as 
a  scattering  problem.  Both  cases  can  be  considered  at  once  by  analyzing 
the  objects  in  an  impressed  field  generated  by  the  source  (J. ,  MJ.  The 
problem  is  depicted  in  Figure  1,  where  the  objects  can  be  arbitrary 

closed  or  open  surface  and  wire  structures,  and  a  wire  may  contact  a  sur¬ 
face. 


Let  S  denote  the  surface  enclosing  the  bodies,  and  n  be  the  unit 
outward  normal  to  S.  In  the  presence  of  the  bodies,  the  sources  (J^ , 
M.)  generate  the  fields  (E,  H).  These  fields  are  considered  to  be 
time  harmonic,  and  the  time  dependence  eJayc  is  suppressed.  From  the 
surface-equivalence  theorem  of  Schelkunoff  [  183,  the  field  interior 
to  the  surface  S  will  vanish  without  changing  the  exterior  field 
if  the  surface  current  densities 

J$  =  n  x  H  (1) 

Ms  =  E  x  n  (2) 

are  introduced  on  the  surface  S,  as  shown  in  Figure  2.  The  null  field 
condition  inside  the  bodies  allows  the  bodies  to  be  replaced  by  the 
ambient  medium  without  altering  the  field  anywhere. 
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SCATTERERS 


Figure  1.  The  scatterers  in  an  impressed  field  generated  by  the  source 

(ii.  M,-). 


Figure  2.  Surface  currents  (J  ,  M  )  are  placed  on  S  without  changing 
the  exterior  fields^  but  making  fields  interior  to  S  vanish. 
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If  ( ,  M.j)  generate  the  incident  fields  (E^,  H,)  in  the  ambient 
medium  (p,e  ),  then  the  scattered  fields  can  be  defined  as 


Hs  =  H  -  H,  (4) 

The  scattered  field  is  radiated  by  the  surface  currents  (Js,  M  ) 
in  the  ambient  medium  (Figure  3).  This  definition  allows  the  total 
fields  (£,  H)  to  be  a  sum  of  the  incident  fields  (E^,  HJ  and  the 
scattered  fields  (£$,  H$). 

A  test  source  (J^,  NL)  is  now  introduced  interior  to  S  to  form  a 
reaction  1  with  the  existing  fields  (£,  H).  Since  the  total  fields 
inside  the  surface  S  are  zero,  the  reaction  of  the  test  sources  and  the 
fields  inside  S  must  be  zero,  i.e. 

(i-Is  -  V«5>  ds  +  #s  -  SU-fli)  1=  -  0-  (5) 

Let  the  fields  radiated  by  (J^,  M^)  in  the  ambient  medium  be  (E^, 

H^).  By  the  reciprocity  theorem,  (5)  can  also  be  written  as 

<JS  -tn  -  «s  -V  *  ♦  US  -I*  -  Mi  •  Hm)  M  *  0  (6) 

V 

where  V  is  the  volume  occupied  by  the  impressed  sources  (J. ,  NT ) . 
Equations  (5)  and  (6)  are  different  forms  of  the  reaction  Integral 
equation  ( R IE )  developed  byRumseyClX  Equation  (6)  was  applied 
by  Richmond  C3  and  5]  for  his  thin  wire  analysis.  The  RIE  is  a  general 
integral  equation  describing  the  electromagnetic  fields.  If  only 
electric  current  test  sources  are  used.  Equation  (6)  reduces  to  the 
widely  known  EFIE  and  the  use  of  only  magnetic  test  sources  will  reduce 
Equation  (6)  into  the  MFIE  [2].  In  this  work,  only  the  electric  test 
sources  are  used.  The  bodies  are  considered  to  be  perfectly  conducting 
and  M$  =  0. 

Although  Equation  (6)  was  derived  with  the  assumption  that  all  sur¬ 
faces  were  closed  so  that  Schekunoff's  surface  equivalance  theorem  [183 
can  be  used,  it  can  also  be  applied  for  open  surfaces.  In  the  case 
of  an  infinitely  thin  flat  plate,  3S  ls  the  vector  sum  of  the  current 
on  the  top  and  bottom  surfaces  if  electric  test  sources  are  used  C203. 


The  following  shows  how  Equation  (6)  is  solved  by  the  moment 
method.  The  unknown  current  _J$  is  expanded  in  terms  of  a  set  of  N  basis 
(expansion)  functions,  Fp: 


is 


N 

y  l  f„  , 

n  =  1  " 


(7) 


and  Equation  (6)  is  enforced  for  N  electric  test  sources  placed  in 
S.  Thus,  Equation  (6)  reduces  to  a  set  of  simultaneous  linear  equations 
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Figure  3.  The  exterior  scattered  fields 


I 


7 


Us.  Hs> 

MEDIUM  (  /x  ,e  ) 

generated  by  (J^,  M$). 


I 


m  =  1,  2  ...  N 


(8) 


1  zmn  "  ^m  ’ 
n=l 

where 

Z__  -  -  f f  E  ‘  F  ds  (9) 

mn  J  J  — m  — n  v  ' 

n 

vm  -  ///  (Ji  •  I*  -  Mi  •  V  dv.  (10) 

V 

A.  L. 

The  integration  in  Equation  (9)  is  over  the  surface  of  the  nLn  expan¬ 
sion  mode. 


8.  Choice  of  Test  Modes 


The  choice  of  the  testing  functions  is  dictated  by  several  factors 
including:  computation  time  for  the  mutual  impedances,  convergence 
rate,  and  the  degree  of  accuracy  required  for  the  final  result.  Obvious 
ly,  the  choice  of  a  delta  function  (point-matching)  over  the  full  sur¬ 
face  test  function  will  significantly  reduce  the  computation  time. 

In  this  work,  the  testing  functions  are  chosen  to  be  the  same  as  the 
expansion  functions  (Galerkin's  method)  for  the  following  reasons: 

a.  Galerkin's  method  will  result  in  a  symmetric  impedance 
matrix.  Thus,  only  half  of  the  impedance  elements, 

i.e.  the  lower  triangular  matrix  (L),  are  calculated  and 
stored.  Also,  solving  for  the  unknowns  (by  the  standard 
Crout  decomposition)  with  the  triangular  matrix  (L)  is 
faster  than  with  the  full  impedance  matrix. 

b.  It  is  felt  that  Galerkin's  method  will  result  in  a  versa¬ 
tile  and  general  computer  code.  Thus,  backscatter,  as 
well  as  numerically  difficult  antenna  problems,  can  be 
solved  by  the  same  general  purpose  computer  code.  In  parti¬ 
cular,  when  considering  the  point-matching  technique  for 

a  wire  attached  to  a  plate,  the  solution  could  be  overly 
sensitive  to  the  proximity  of  the  attachment  point  and 
the  match  point.  However,  experience  with  rectangular 
surface  patch  modes C 16 □  indicated  that,  with  surface  test¬ 
ing,  the  attachment  point  may  be  anywhere  in  the  surface 
patch  grid. 

c.  Galerkin's  method  will  ordinarily  converge  faster  than  point 
matching.  Thus,  for  a  given  surface  patch  grid,  the  result 
obtained  by  Galerkin's  method  will  be  closer  to  convergence 


8 


than  with  point  matching.  This  also  implies  that  a  smaller 
number  of  unknowns  is  required  for  a  given  problem. 


C.  Expansion/Testing  Functions 

In  order  to  model  a  geometry  with  thin-wires,  flat  surfaces, 
and  wire-surface  connections,  three  basic  types  of  modes  (wire  dipole 
modes,  surface  dipole  modes,  and  the  wire-surface  attachment  modes] 
are  used.  Singly  curved  surfaces  can  also  be  modelled  using  piecewise 
flat  approximation.  All  of  the  modes  will  have  a  sinusoidal  variation 
with  free  space  wave  number,  k,  in  the  direction  of  current  flow. 


1.  Thin-wire  Mode 

The  wire  mode  used  is  the  piecewise-sinusoidal  (PWS)  V-dipole 
consisting  of  two  sinusoidal  monopolos  developed  for  the  thin-wire  formu¬ 
lation  by  Richmond C  19] .  Figure  4  shows  a  V-dipole  with  a  180°  internal 
angle  lying  on  the  z-axis.  Using  the  thin-wire  approximation,  the 
current  on  this  dipole  is  given  by 


iw  -  2 

is  '  TYa 

where  P^  and  P^  represent  pulse  functions  with  unit  value  when  z^<  z<  z2 
and  Zp<z<  z^,  respecti vely,  and  are  zero  elsewhere.  Also,  a  is  the 
wire  radius  and  k  =  2"  /A  .  These  modes  are  placed  in  an  overlapping 
array  on  the  wire  to  ensure  continuity  of  current  on  the  wire.  The 
advantage  of  PWS  mode  is  that  the  near-zone  fields  and  the  mutual  imped¬ 
ance,  zmn,  are  known  in  closed  form  [  22]. 


sin  k  (z-z^)  sin  k  (z-^-z) 

P1  sin  k  (zrzj7  +  P2  sin  k  (zyZp) 


(11) 


2.  Surface-patch  Mode 

The  general  quadrilateral  surface-patch  mode  is  a  surface  V-dipole 
consisting  of  two  quadrilateral  piecewise-sinusoidal  surface  monopoles. 

A  surface  V-dipole  with  interior  angle  of  180°  is  shown  in  Figure  5. 
Consider  some  general  point  in  monopole  A  (shown  as  the  heavy  dot  in 
Figure  5),  which  has  terminal  length  V  and  end  length  U.  Construct 
a  line  through  the  point  which  intersects  the  end  and  terminal  such 
that  u/v  =  U/V.  Let  L  be  the  length  of  this  line  and  i  be  the  distance 
along  the  line  from  the  end  to  the  point.  Then  the  surface  current 
density  on  monopole  A  is 


Ja  = 


C  sin  kH  : 

~wTr/LT  1 


(12) 


9 


MONOPOLE  A 


where  W( H  /L)  is  the  local  width  of  monopole  A,  and  C  is  a  normalizing 
constant  such  that  the  terminal  current  is  one  ampere.  The  surface 
current  density  on  monopole  B,  Jn,  is  defined  in  a  similar  way,  but  with 
the  current  flowing  in  the  -j i  direction,  this  mode  is  a  general ization 
of  the  rectanqu'iar  PWS  surface-patch  dipole  used  by  Wang  £133  as  well  as 
Newman  and  Pozar[16]  . 

3.  Attachment  Mode 

The  attachment  mode  is  used  to  establish  continuity  of  current 
at  the  wire/surface  junction  and  to  ensure  that,  in  the  immediate  vicin¬ 
ity  of  the  attachment,  the  surface  current  density  has  the  proper  p 
polarization  and  1/p  dependence.  This  mode  consists  of  two  parts: 
a  wire  monopole  and  a  disk  monopole,  as  shown  in  Figure  6.  The  wire 
monopole  current  density  is  similar  to  the  thin  wire  monopole  mode: 


^s 


I  sin  k (z^-z) 


2tt  a  sin  kz. 


z. 


(13) 


The  disk  current  density  is 

iD  _  -sin  k(b-P  ) 

— s  2 up  sfn  k(b-aT 


P  ,  a  <  P  <  b, 


(14) 


where  a,b  are  the  inner  and  outer  radii  of  the  annulus,  respectively. 
Note  that  the  total  current  on  the  disk  at  p-a  >s  equal  to  the  total  cur 
rent  on  the  wire  segment  at  z=0  insuring  continuity  of  current  at  the  at 
tachment.  Also,  observe  that  the  sin  k(b-p)  function  in  the  numerator 
of  Equation  (14)  forces  the  disk  current  to  be  zero  at  the  disk  edge, 
p  =b.  It  is  this  property  which  allows  the  disk  to  be  placed  on  the 
surface  and  still  maintain  continuity  of  current  on  the  surface.  The 
sin  k(b-  p)  function  was  chosen  (rather  than,  say,  b-p  or  sin  (b-  p) ) 
since  this  permits  the  fields  of  the  disk  to  be  obtained  with  only 
one  numerical  integration. 

The  attachment  mode  is  used  by  placing  it  directly  over  the  sur¬ 
face-patch  modes  where  the  wire  meets  the  surface.  There  is  no  special 
way  in  which  the  disk  must  be  situated  with  respect  to  the  surface- 
patch  modes;  that  is,  the  disk  may  be  centered  at  midpoints,  corners, 
or  other  points  on  the  surface-patch  modes.  The  only  restriction  is 
that  the  attachment  point  be  at  least  0.1  X  from  the  edges  of  plates. 

The  outer  radius  of  the  disk  mode,  b,  should  be  in  the  range  from  about 
0.1  X to  0.21AC163. 

For  antennas,  feeds  may  be  placed  at  the  endpoints  of  any  wire 
segment  including  the  attachment  segment.  For  example,  the  feed  of 
a  monopole  on  a  finite  ground  plane  may  be  modelled  by  inserting  a 
generator  at  the  base  of  the  monopole  between  the  attachment  mode 


Figure  6 
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Attachment  dipole. 
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wire  and  disk  monopole.  Delta-gap  reed  is  used  in  this  work,  although 
a  magnetic  frill  model  could  also  be  used. 


4.  Plate  Overlap  Mode 

The  overlap  modes  allow  a  non-zero  continuous  current  at  the 
plate-to-plate  junction.  They  are  identical  to  the  surf ace-patch  dipole 
mode  discussed  earlier.  The  edges  .if  the  overlap  surface  dipoles  need 
not  coincide  with  the  edges  of  the  surface  dipole  modes  on  either  plate, 
thus  allowing  the  intersection  of  plates  of  different  sizes. 


D .  Computation  of  Impedance  E  i enients 

Many  computationally  efficient  techniques  developed  for  the  evalu¬ 
ation  of  the  mutual  impedances  associated  with  rectangular  surface- 
patches  L 21 3  can  also  be  used  for  the  quadr i lateral  surface-patches. 

This  section  will  discuss  these  techniques  and  how  they  are  used  in 
finding  the  mutual  impedances  of  quadri lateral  surf ace-patches.  The 
majority  of  the  mutual  impedances  are  either  the  surface-to-surface 
dipole  or  disk-to-surface  dipole  impedances.  The  disk-to-surface  dipole 
impedance  is  calculated  in  the  same  way  as  described  in  ref  erenceC  213. 
The  surface-to-surface  dipole  impedance  is  discussed  below. 

The  general  expression  for  the  impedance  matrix  element  given 
in  Equation  (9)  is  evaluated  numerically  by  considering  the  surface- 
dipole  modes  as  made  up  of  PWS  filaments.  The  surf ace-to-surf ace  dipole 
impedance  is  a  sum  of  four  surface-patch  monopole-to-monopole  impedances, 

i.e. 


mn 


ztAeA  +  ztAeB 


■f 


t'3eA  t3e8 
z  +  z 


(15) 


where  the  superscript  tA  and  e^  denote  the  test  monopole  A  and  expansion 
monopole  8,  respectively.  Also,  to  and  e„  are  defined  similarly. 

Each  surface-patch  monopole-to-monopole  impedance  is  evaluated  by  sum¬ 
ming  up  all  the  fi lament-to-fi lament  impedances C  213.  If  the  test 
and  expansion  surface  monopoles  are  each  represented  by  M  filaments, 
then,  in  general,  each  monopole-to-monopole  impedance  requires  M^  evalua¬ 
tions  of  the  fi lament- lo-fi  lament  impedances.  If  both  monopoles  are 
rectangular,  and  if  either  the  surface  current  directions  are  parallel 
or  the  vector  transverse  to  the  surface  monopoles  are  parallel,  then 
only  2M  of  these  impedances  are  evaluated  [  213. 

Most  of  the  time  spent  in  filling  out  the  impedance  matrix  is 
used  in  finding  the  fi lament-to-fi lament  impedances.  Thus,  computation 
time  depends  on  the  number  of  filaments  used  to  represent  a  surface 
monopole.  In  evaluating  the  self  impedance  or  mutual  impedance  of 
two  closely  spaced  surface  monopoles,  experience  has  shown  that  a  fila¬ 
ment  spacing  of  0.07  k  or  about  five  filaments  per  0.25k  wide  surface 
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patch  is  adequate  for  representing  the  surface  patch.  More  numerically 
difficult  problems,  such  as  finding  near  field  quantities  of  an  antenna, 
may  require  a  slightly  larger  number  of  filaments  usually  no  more  than 
seven  filaments  per  0.25X  wide  surface  patch.  Simpson's  rule  of  inte¬ 
gration  is  used  to  sum  up  the  fi lament-to-fi lament  impedances  for  a 
surface  monopole-to-monopole  impedance.  Usually,  when  the  monopoles  do 
not  intersect  each  other,  a  three-point  Simpson's  rule  (or  repre¬ 
senting  the  surface  patch  by  three  filaments)  is  used.  For  surface 
monopoles  spaced  more  than  0.25A  apart,  a  one-filament  representation 
for  each  surface  monopole  is  used.  As  an  example,  the  reactance  of 
two  rectangular  monopoles  is  plotted  versus  the  separation  distance, 

0,  in  Figure  7.  Plotted  are  the  exact  reactance  (solid  line)  and  the 
reactance  obtained  by  representing  the  surface  monopole  by  one  filament 
on  the  center  of  the  monopole  (dotted  line).  The  difference  of  the 
two  curves  results  from  the  phase  and  amplitude  variations  of  the  fila- 
ment-to-fi lament  impedance  as  the  integration  is  done  across  the  surface 
patct|gjj.  These  phase  and  amplitude  variations  can  be  approximated  by 

^  where  R  is  the  average  distance  of  the  two  filaments.  To 

obtai^a  more  accurate  result  by  the  one-filament  representation, 

— )  is  used  in  the  following  manner.  If  is  the  filament- 

to-fi lament  impedance,  then  the  variation  due  to  the  nth  expansion 
surface-patch  monopole  is  taken  into  account  by: 


cfn  =  hi 


Re 


jkR 


-Mi 


4e 


-jkR 


-jkR, 


(16) 


where  R1,  R,  and  R?  are  the  average  distances  from  the  test  filament 
to  the  closer  edge,  center,  and  further  edges  of .the  surface-patch 
monopole.  Similarly,  the  variation  due  to  the  nr  test  surface-patch 
monopole  is  accounted  for  by: 


mn 


!flL  Re 


jkR 


(17) 


The  imaginary  part  of  zmn  is  plotted  as  dashed  lines  in  Figure  7  and 
is  very  close  to  the  accurate  impedance  value  (solid  line).  The  real 
part  of  zmp  also  behaves  in  the  same  manner. 

Another  time-saving  technique  is  the  use  of  an  equivalent  reactance 
when  two  PWS  filaments  are  very  close  to  each  other.  This  is  used 
only  when  two  surface-patches  coincide  or  share  a  common  side.  The 
equivalent  reactance  is  derived  in  the  same  manner  as  is  done  in  refer¬ 
ence  C 213.  The  reactance  of  two  closely  spaced  filaments  can  be  written 
as 


Y(x)  =  S.n  (x) 


(18) 
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REACTANCE  (OHMS) 


The  reactance  of  two  rectangular  PWS  surface-patch  monopoles 
M  denotes  the  number  of  PWS  filaments  used  to  represent  the 
surface  patch. 


Jr 


where  Y  is  the  reactance  and  x  is  the  separation  distance.  C,  and  C~ 
are  constants  which  can  be  evaluated  by  c 

Y( A  x)  -  Y (—7^) 

r*  _  ^  /  1  r\  \ 


C,  =Y(Ax)  -  Co  £  r>  (Ax) 


Equation  (18)  can  be  integrated  analytically  to  give  an  equivalent 
reactance  at  x=0: 


Y(0)  =  -Y(A  x)  +  2Cj  +  2C2  (£  n  (  A  x)  -  1) 


Although  Equations  (18)  through  (20)  were  derived  for  parallel  filaments, 
it  is  also  applicable  to  nonparallel  filaments  as  shown  in  Figure  8. 

Using  the  equivalent  reactance  at  x=0  (Equation  21)  avoids  the  £ n- 
singularity  and,  thus,  signif icantly  increases  the  convergence  rate 

12  Q. 

Another  factor  to  keep  in  mind  when  finding  the  PWS  filament- 
to-fi lament  impedance  is  whether  to  evaluate  the  closed  form  expression 
[221]  or  to  use  the  less  time  consuming  numerical  integration.  This 
problem  is  addressed  in  Appendix  A. 
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CHAPTER  III 


AUTOMATIC  GENERATION  OF  SURFACE  PATCH  MODE  LAYOUT 


The  problem  of  setting  up  surface-patch  modes  over  an  arbitrary  pol¬ 
ygonal  flat  plate  is  considered  in  this  chapter.  Three  different 
subdivision  schemes  are  discussed.  The  first  scheme  is  the  Gordon- 
Hall  (G-H)  procedure  L 23] ,  and  the  last  two  are  procedures  devised 
to  achieve  an  optimal  mode  set-up  by  satisfying  the  criteria  given 
in  the  next  section.  Examples  are  used  to  illustrate  the  appl  icabi  1  i  ty 
of  each  scheme. 


A.  Introduction 


The  set  of  coordinates  describing  all  the  surface  dipole  modes 
forms  a  large  portion  of  the  input  data  required  for  the  moment  method. 
Although  coordinates  for  each  mode  can  be  manually  entered  into  a 
moment  method  computer  code,  it  is  desirable  that  the  mode  set  up 
be  done  automatically  with  the  input  being  the  geometry  of  the  plates. 

A  reasonable  mode  generation  algorithm  should  minimize  the  user  data 
input  and  be  consistent  with  the  following  criteria: 

a.  In  general,  there  must  be  two  orthogonal  sets  of  surface 
dipole  modes  and  their  linear  combination  can  approximate 
the  actual  plate  current. 

b.  At  the  edge  of  the  plate,  the  component  of  current  perpend¬ 
icular  to  the  edge  must  be  zero  and  that  parallel  to  the 
edge  must  be  non-zero. 

c.  The  sides  of  modes  at  the  boundary  of  the  plate  must  match 
the  curvature  of  the  boundary  adequately,  especially  on 
sharp  bends. 

d.  The  size  of  most  modes  should  be  about  the  same  as  a  speci¬ 
fied  nominal  size  and  the  length  of  each  monopole  should 
not  exceed  some  maximum  size. 

e.  Where  possible,  the  shape  of  the  monopoles  should  be  rectang¬ 
ular  and  the  directions  of  currents  in  the  majority  of  the 
modes  should  be  parallel. 


Criterion  (a)  is  required  for  the  solution  to  converge,  while 
criterion  (b)  is  a  relaxed  requirement  of  the  edge  condition  C 24 J. 

This  condition  states  that,  in  the  vicinity  close  to  the  edge  of  an 
infinitely  thin  plate,  the  currents  perpgndiculaf  and  parallel  to 
the  edge  have  the  asymptotic  behaviors  p  2  and  p  2,  respectively,  where 
Pis  the  perpendicular  distance  from  the  edge.  However,  the  calcula¬ 
tions  of  the  input  impedance  of  a  monopole  on  a  square  plate  by  Pozar 
and  Newman[25j  suggested  that  the  surface  patch  dipoles  adequately 
approximate  this  behavior  a  least-mean-square  sense,  and  special  edge 
modes  with  the  required  asymptotic  behaviors  do  not  s ignif icantly 
improve  the  final  solution.  On  the  contrary,  k ichmond  E 263  has  added 
the  edge  modes  (with  proper  behavior  at  the  edge)  to  the  other  basis 
functions  for  strip  grating  and  found  that  the  convergence  is  improved. 
For  simpl icity, these  special  edge  modes  have  been  avoided,  and  satisfact¬ 
ory  results  (see  Chapter  Four)  have  been  obtained.  Criterion  (c) 
ensures  that  the  boundary  condition  is  satisfied  on  the  boundary  of  the 
plate. 

Although  criteria  (d)  and  (e)  are  not  necessary  requirements, 
they  are  very  important  for  the  best  computational  efficiency.  By  keep¬ 
ing  the  sizes  of  the  patches  as  close  to  the  nominal  value  as  possible, 
as  stated  in  criterion  (d),  the  total  number  of  unknowns  can  be  minimiz¬ 
ed.  Rectangular  modes  are  desirable  by  virtue  of  their  rapid  converg¬ 
ence  rate  and  fast  self  and  mutual  impedance  evaluation C 213.  The 
expression  for  the  mutual  impedance  between  two  closely  spaced  parallel 
current  filaments  is  simpler  to  evaluate  than  that  of  the  non-parallel 
case;  thus,  parallel  filaments  result  in  reduced  overall  computation 
time. 


The  corners  (or  nodes)  and  sides  of  the  surface  dipole  modes 
form  a  grid  or  mesh  which  can  be  automatically  generated.  Okon  and 
Harrington C 27 3 1 isted  several  works  done  on  automatic  mesh  generation 
and  described  the  Gordon-Hall  procedure,  which  is  relatively  easy 
to  implement  and  satisfies  several  of  the  criteria  listed  previously. 


B.  Gordon-Hall  (G-H)  Procedure  C233 


The  G-H  procedure  forms  a  curvilinear  coordinate  system  in  an 
arbitrary  region,  R,  by  a  one-to-one  mapping  (U:  S  -*•  R)  from  a  unit 
square  A'B'C'O'  in  the  s-t  plane  (S  domain)  to  the  region  bounded 
by  ABCD  in  the  x-y  plane  (R  domain)  (see  Figure  9).  In  the  following 
discussion,  R  is  a  planar  surface,  but  the  method  is  applicable  to 
an  arbitrary  curved  surface  embedded  in  three-dimensional  space. 

The  theory  assumes  the  existence  of  a  continuous  vector-valued  function 
F  which  maps  S  onto  R.  However,  in  practice,  the  geometric  description 
of  R  is  given  in  terms  of  its  boundary  and  F  becomes 


F  (s,t) 


x  (s,t)' 
y  (s,t) 


(22) 


where  x  (s,t)  and  y  (s,t)  are  parametric  equations  describing  the 
boundary  of  R  (3R)  in  the  x-y  plane.  A  class  of  mapping  or  function 
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U  which  interpolate  £  at  a  non-denuwerablo  set  of  points  (transfinite 
property)  comprising  the  boundary  oi  S  is  given  by  Gordon  and  Hall 
f  23 1  as: 

U  ( s ,  t )  =  ( 1-s  )F(0,  t)  +  s  T(l,t)  +  ( 1-t.  )F(s,0)  i  t  F  ( s ,  1 )  - 

( 1-s)  ( 1- 1 )F (0,0)  -  ( 1  - s ) t.F (0,1)  -  s(l-t)F(l,0)  - 

stF(I.I)  (23) 

U  can  be  described  as  a  Function  which  induces  generalized  coodinates 
s  and  t  onto  the  region  K.  ihis  i  an  best  be  illustrated  by  examples 
given  below. 


1.  The  Quadri  lateral  [;'/.] 


The  steps  in  forming  a  d-ii  aiesh  generator  tor  a  quadrilateral 
are  ( s>  e  Figure  10) : 


a.  Choose  Four  points  on  DF  which  will  separate  the  boundary 
of  the  quadr :  lafni  ,j  i  in,.;  four  segments.  Tliese  segments 
should  he  approximate ly  equal  in  order  to  have  a  uniformly 
spaced  mesh.  In  this  example,  it  is  obvious  that  the  four 
points  are  the  corners  of  tii^  quadrilateral. 


b.  The  four  segment >  obtained  above  are  to  be  the  images  of 
the  four  sides  of  the  unit  square  S  under  the  mapping  F. 
That  is,  the  segments  are  graphs  of  the  four  vector-valued 
functions:  for  6  5  s,t  C  1, 


F(s,0) 

F(l,t) 

F(s,l) 

F(0,t) 


s(k2-k1)\^ 

s(y?-y,  )J 

t(*3-x?)\ 

t(y3-.y?)/ 

s( x3"x4 

5(y3-.y,/|)/ 

t(x4"x])> 
t  ( v/, -y  i  >/ 


(24) 

(25) 

(  2b) 

(2  7) 


c. 


The  transfinite  map  H  (or  the  b-li  mesh  generator)  is  construct¬ 
ed  by  substituting  equations  (24)  through  (27)  into  Equation 
(p 3),  which  result/,  in: 


r 

i 


U  (s,t) 


/x(s,t)\ 
\y(s,t)  J 


/V'l\  7TW4\  t 

\y2‘yl  /  \yry2+y3-y4/ 


t 


(28) 


A  grid  formed  by  (28)  is  shown  in  Figure  11. 

A  special  case  of  a  quadri lateral  is  the  triangle.  The  fourth 
corner  is  arbitrarily  set  at  the  midpoint  of  its  longest  side.  The 
mesh  formed  by  Equation  (28)  is  shown  in  Figure  12. 


2.  The  El  1  ipse  C  J 

This  example  differs  from  the  previous  ouo  in  two  '’aspects:  first, 
it  does  not  have  a  set  of  four  well  defined  points  for  segmenting 
aR  and,  second,  aR  is  a  curved  line  instead  of  straight  line  segments 
(see  Figure  IF).  However,  tiie  basic  steps  in  finding  U  are  the  same: 


a.  The  important  consideration  in  choosing  the  four  corner 
points  is  to  choose  them  so  that  the  resulting  arc  lengths 
are  about  the  same.  Thus,  the  four  points  (a,0),  (o,b), 
(-a,0),  and  (0,-b)  are  chosen. 

b.  The  function  F  is:  for  0  <s,t  <1, 


F(s,0) 


( 


S  IT 

6  cos  — 2 — 


b  sin 


S  71 

~T~ 


(29) 


F(l,t) 


a  cos  4  (t+1) 
b  s in  4  ( t  +  l ) 


(30) 


F ( s,  1 ) 


F(0,t) 


(a  cos 

-b  sin 


(FI) 


(32) 
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Figure  11.  The  curvilinear  coordinates  formed  by  the  G-H  Procedure  for 
the  quadri lateral  in  Figure  10. 


2<\ 


c. 


Substituting  Equations  (29)  through  (32)  into  (23)  yields: 


U(s,t)  = 


^(l-s)cos^  -s  sinr-|+(l-t)cos^-|  -t  sirv|+s+t 
[s  cos^|-(l-s)sirA|+(l-t)sin^  -t  cos^|+t-s 


(33) 


Note  that,  in  this  case,  U  maps  the  straight  lines  in  the 
s-t  plane  into  curved  lines  in  the  x-y  plane,  as  shown  in 
Figure  14.  Since  only  straight  line  segments  are  allowed 
for  the  surface  patch  mode,  the  mesh  is  formed  by  connecting 
the  intersection  points  by  straight  line  segments  as  shown  in 
Figure  15. 

G-H  procedure  is  an  automatic  mesh  generator  which  satisfies 
cr-teria  (a)  through  (c)  in  section  A  of  this  chapter.  However ,  criter¬ 
ia  (d)  and  (e)  are  usually  not  met  for  an  arbitrary  boundary.  In 
order  to  have  criterion  (d)  satisfied,  the  four  arcs  describing  the 
3R  must  be  carefully  chosen  such  that  they  are  approximately 
of  equal  length.  Note  that  each  arc  can  be  a  straight  line,  curved 
lines,  or  two  or  more  straight  line  segments  f?3”.  However,  a  bad 
choice  of  the  arc  sections  may  result  in  a  mapping  U  which  does  not 
map  S  onto  the  correct  planar  surface  of  R  r2 .?•!.  To  alleviate  this 
problem,  a  region  R  with  complex  boundary  may  be  subdivided  into  subre¬ 
gions  and  each  region  is  mapped  independently.  This  will  require 
the  user  to  parameterize  the  boundaries  of  several  subregions,  and  these 
subregions  must  be  joined  electrically  by  overlapping  surface-patch 
dipole  modes.  Since  criterion  (e)  is  seldom  met  in  an  arbitrary  region 
R,  this  results  in  inefficient  evaluation  of  the  impedance  which  trans¬ 
lates  into  a  large  amount  of  computation  time.  Alternative  subdivision 
schemes  are  proposed  to  find  a  more  desirable  mesh  for  regions  R  which 
are  polygons. 


C.  Subdivision  Method  1 


The  objective  of  this  section  and  the  next  is  to  present  algorithms 
such  that,  when  implemented  in  the  form  of  a  computer  code,  can 
provide  the  coordinates  for  a  set  of  surface  dipole  modes  which  satis¬ 
fies  all  the  criteria  given  in  Section  A.  it  is  assumed  that  an 
arbitrarily  shaped  plat j  can  be  adequately  described  as  a  polygon. 

Thus,  the  plate  can  be  described  by  a  set  of  points,  ordered  in  either 
clockwise  or  counter-clockwise  direction,  and  the  straight  lines  joining 
these  points  form  the  boundary  of  the  plate.  The  basic  approaches  used 
in  these  two  subdivision  schemes  are  discussed  in  Appendix  C. 

The  general  approach  of  this  scheme  is  to  choose  a  side  and  divide 
the  plate  into  strips  which  are  parallel  to  this  side.  The  width  of 
the  strips  is  to  be  about  the  same  or  less  than  a  specified  maximum 
size,  W.  Each  end  side  of  the  strips  should  be  a  straight  line  segment 
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gure  14.  The  curvilinear  coordinate  system  induced  by  the  transf 
ite  mapping  U. 


gure  15.  A  mesh  for  the  overlapping  surface-patch  dipole  modes. 


which  is  a  part  of  the  plate  boundary.  A  grid  can  be  formed  by  segment¬ 
ing  the  strips  into  an  equal  number  of  blocks  and  the  maximum  length 
for  the  blocks  should  be  less  than  or  equal  to  the  maximum  width, 

W.  The  details  of  the  algorithm  is  given  below  (see  Figures  16  and 
18  -  note  that  the  arrows  indicate  the  direction  of  the  overlapping 
surface  patch  dipole  modes): 


a.  Choose  a  side  such  that  strips  parallel  to  this  side  (refer¬ 
ence  side)  can  be  formed  throughout  the  plate.  Usually 

the  longest  side  is  chosen.  As  shown  in  Figure  16,  either 
sides  A8,  DE,  or  AE  can  be  chosen,  but  choosing  other  sides 
will  not  allow  the  plate  to  be  divided  into  parallel  strips. 
That  is,  the  reference  side  must  be  chosen  such  that  the 
perpendicul ar  vectors  from  the  reference  side  to  all  other 
corners  of  the  plate  must  all  point  in  the  same  direction. 

Also,  three  consecutive  vectors  should  have  the  property  that 
the  length  of  the  middle  vector  is  not  the  shortest  among  the 
three  vectors.  That  is  the  situation  as  shown  in  Figure  17, 
where  AE  chosen  as  the  reference  side  is  not  permissible.  If  AB 
is  not  parallel  to  DE,  then  the  strip  with  a  side  coinciding 
to  DE  will  not  be  trapezoidal.  The  strip  width  is  found 
by  dividing  the  perpendicular  distance  from  the  reference 
line  (AB  in  this  example)  to  the  nearest  corner  by  the  small¬ 
est  integer  ft  that  will  make  the  result  less  than  or  equal 
to  W.  Note  that  strip  widths  can  vary  depending  on  the 
distance  from  the  reference  side  to  the  corners.  If  the 
length  DE  is  zero,  then  the  strip  with  side  DE  will  have 
very  short  and  zero  length  current  elements.  This  is  avoided 
by  reducing  the  width  of  this  particular  strip  so  that  there 
will  be  no  side  with  zero  length.  Since  there  will  be  current 
filaments  along  the  end  sides  of  each  strip  (see  Figures 
16  and  18),  each  strip  must  have  no  end  sides  whose  lengths 
exceed  0.25  X  If  this  is  not  satisfied,  the  strip  must  be 
further  divided  into  smaller  strips. 

b.  Find  the  minimum  integer,  fL,  which  the  longest  strip  can 
be  segmented  such  that  the  length  of  each  section  along 
this  strip  does  not  exceed  W.  Each  strip  is  then  divided 
into  sections,  resulting  in  a  mesh  as  shown  in  Figures 
16  and  18.  Note  that,  if  the  two  end  sides  of  a  strip  are 
not  parallel,  only  trapezoidal  blocks  can  be  formed. 

The  mesh  formed  by  this  scheme  will  have  about  half  of  the  surface 
dipoles  with  current  filaments  that  are  parallel  to  each  other  (see 
Figures  16  and  18).  This  is  due  to  the  fact  that  usually  only  trape¬ 
zoidal  blocks  can  be  formed.  However,  the  more  disturbing  feature 
of  this  scheme  is  that,  if  there  are  strips  whose  lengths  are  much 
shorter  than  others,  there  will  be  a  relatively  large  number  of  very 
small  (or  short)  surface  dipoles  in  the  shorter  strips.  In  spite 
of  these  shortcomings,  this  scheme  is  still  attractive  is  a  plate 
can  be  divided  into  strips  whose  lengths  do  not  vary  significantly. 
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gure  16.  A  set  of  overlapping  surface-patch  dipole  modes  formed 
side  a  fin-shaped  polygon  by  using  subdivision  method  1 
The  arrows  indicate  the  direction  of  the  currents. 
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Figure  18.  A  set  of  overlapping  suri 
side  a  regular  octagon  b' 
area  inside  the  polygon 
a  0.3  X  radius. 


A  modification  of  this  algorithm  results  in  the  next  scheme  where 
two  meshes  are  formed  for  the  two  current  polarizations. 

D.  Subdivision  Method  2 

This  subdivision  scheme  is  similar  to  the  previous  one  in  that 
strips  are  formed  from  one  side  of  the  plate.  However,  the  strips 
are  not  segmented  by  the  same  number  of  sections  N t ;  instead,  each 
strip  is  segmented  into  a  minimum  number  of  sections.  This  allows 
the  sizes  of  the  surface  monopoles  to  be  more  uniform.  The  interior 
sides  are  adjusted  such  that  most  of  the  segments  in  a  strip  will 
be  rectangular  (see  Figures  21  and  22).  This  choice  of  segmentation 
will  allow  only  one  set  of  surface  dipole  currents  in  the  direction 
along  the  length  of  the  parallel  strips.  The  other  set  of  dipole 
modes  are  formed  by  dividing  the  plate  into  strips  again,  but  this 
time  the  strips  are  perpendicular  to  the  previous  set  of  strips  (see 
Figure  23).  The  details  of  this  subdivision  scheme  are: 

a.  This  step  is  similar  to  step  (a)  in  subdivision  method  1. 
However,  relatively  short  strips  are  allowed.  It  is  possible 
that  a  plate  geometry  (such  as  that  shown  in  Figure  17) 

will  not  allow  only  one  strip  with  the  ends  being  the  perimet¬ 
er  of  the  plate  to  be  formed.  To  get  around  this  problem, 
another  side  may  be  chosen  for  the  reference  side.  Alternat¬ 
ively,  the  plate  can  be  divided  into  two  plates  and  surface- 
patch  dipoles  set  up  for  each  plate.  These  two  plates  are 
then  joined  by  the  plate  overlap  modes. 

b.  Each  strip  is  segmented  into  the  minimal  number  of  sections 
such  that  the  length  of  each  section  should  not  exceed  the 
specified  nominal  value,  W.  If  the  two  end  sides  are  parallel, 
i.e.  AD | ( BC  in  Figure  19,  the  sections  will  be  rectangular. 

If  AD  is  not  parallel  to  BC  and  the  length  of  the  longest 
side  of  the  strip  is  less  than  three  times  the  maximum  width, 

W,  the  strip  is  divided  as  shown  in  Figure  20. 

In  the  case  where  AD  is  not  parallel  to  BC  and  the  length  of 
the  longest  side  is  greater  than  three  times  the  nominal 
length,  an  attempt  is  made  to  segment  the  strip  such  that  the 
segments  in  t hr  middle  of  the  strips  are  rectangular.  As 
shown  in  Figure  21,  if  |AD|<  W  ^/2 ,  then  the  strip  can  first 
be  segmented  by  the  line  A'D'.  '’The  line  A'D'  is  placed  such 
that  the  length  of  AA‘  is  the  maximum  width,  W.  If  the  length 
of  DO'  is  less  than  0.2  W,  the  point  E',  which  is  the  midpoint 
of  segment  DE,  is  used  to  segment  the  strip.  This  will  allow 
the  quadri lateral  AA'E'D  to  have  a  more  uniform  side  length. 

The  same  is  done  for  the  other  end  and  the  segmentation  of  the 
rectangle  EFB'C'  will  result  in  rectangular  sections,  as  shown 
in  Figure  21.  If  } AD )  >  W  ^ j  2  (see  Figure 
22),  then  the  strip  can  first  be  segmented 
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Figure  19.  Subdivision  of  a  rectangular  strip. 


D  C 
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Figure  20.  Subdivision  of  a  trapezoidal  strip. 


Figure  21.  Subdivision  of  a  trapezoidal  strip  with  rectangular  sec¬ 
tions  in  the  middle  part  of  the  strip. 


Figure  22.  Subdivision  of  a  trapezoidal  strip  with  very  long  end  sec¬ 
tions  (AD  and  AD' ) . 
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by  the  line  A'D.  The  line  A'D  is  the  perpendicul ar  line 
from  point  D  to  the  line  segment  AB.  The  line  AD  is  segment¬ 
ed  into  N  segments  such  that  no  section  is  greater  than 
W.  Note  that  line  AA‘  is  segmented  one  segment  less  than 
the  line  AD.  This  allows  the  quadr i lateral  APQD'  to  be 
formed  at  the  end  of  the  strip.  The  line'  segment  B * C *  is 
found  like  the  above  and  the  rectangle  A  B'C'D  is  segmented 
into  rectangular  segments. 

c.  To  form  the  other  set  of  currents  orthogonal  to  the  first 
set,  a  line  is  drawn  from  the  furthest  corner  and  perpendicu¬ 
lar  to  the  reference  side.  Strips  are  then  formed  parallel 
to  this  line  and  the  division  is  the  same  as  that  of  steps 
(a)  and  (b). 

d.  As  shown  in  Figure  23,  there  may  be  some  sides  of  an  arbitrar¬ 
ily  shaped  plate  which  do  not  have  non- zero  currents  parallel 
to  the  edges,  fxtra  strips  are  formed  parallel'  to  these 
edges.  I  he  widths  of  these  strips  are  the  average  width 

of  the  previously  formed  strips  or  smaller. 


The  three  subdivision  schemes  discussed  in  this  chapter  are  just 
a  few  of  the  possible  mode  layouts.  Experience  has  shown  that  subdivi¬ 
sion  methods  1  and  2  are  applicable  to  a  wide  range  of  plate  shapes. 

A  computer  subroutine  has  been  written  which  implements  method  1  and 
method  2.  The  subroutine  chooses  method  1  if  the  difference  of  the 
longest  and  shortest  strip  formed  is  less  than  4W,  and  chooses  method 
2  otherwise.  Also,  if  the  length  of  the  shortest  strip  is  less  than 
0.1  W,  then  method  2  will  be  chosen.  The  user  specifies  the  corners 
of  the  polygon  and  the  maximum  segment  size,  W.  The  subroutine  gener¬ 
ates  the  coordinates  of  the  surf ace-patch  dipole  modes  on  the  polygon 
plate.  The  accuracy  and  efficiency  of  these  mode  set-ups  in  the  moment 
method  is  discussed  in  the  next  chapter  on  numerical  results. 
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CHAPTER  IV 


NUMERICAL  RESULTS 


In  this  chapter,  data  for  the  input  impedance  of  a  monopole  anten¬ 
na  centered  on  a  disk  and  scattering  data  from  different  polygonal 
plates  are  presented.  The  surface-patch  modes  are  formed  automatically 
by  the  subdivision  algorithms  given  in  the  preceding  chapter.  A  maxi¬ 
mum  segment  size  W=0.25A  is  used  for  all  plate  subdivisions.  The  re¬ 
sults  are  compared  with  measurements  and  independent  calculations. 

All  the  computations  were  done  on  the  Digital  Equipment  Corporation 
VAX  11/780,  which  is  three  to  four  times  slower  than  the  IBM-370/168. 


A.  Monopole  Antenna  on  a  Disk 

The  study  of  the  monopole  antenna  on  a  circular  disk  is  of  inter¬ 
est  since  it  leads  to  the  understanding  of  the  effect  of  the  finite 
ground  plane  on  antenna  input  impedance  and  radiation  pattern.  Richmond 
L283  listed  the  work  on  this  subject  and  presented  a  moment  method 
solution.  PWS  current  forms  were  used  for  the  wire  and  the  disk  dipole 
mode  consisted  of  two  concentric  annular  zones  with  radially  directed 
currents  across  the  two  zones. 

In  order  to  use  the  techniques  developed  in  this  work,  the  circular 
disk  is  first  approximated  by  a  polygon.  Thin-wire,  attachment,  and 
surface-patch  modes  are  used  to  represent  the  surface  currents.  The 
surface-patch  modes  are  formed  by  using  subdivision  method  1.  The 
surface  dipole  modes  for  an  octagonal  disk  with  a  radius  of  0.3A  are 
shown  in  Figure  18.  The  input  resistance  and  reactance  of  a  A/4  mono- 
pole  versus  disk  radius  are  shown  in  Figures  24  and  25,  respectively. 
Data  are  shown  for  square  and  octagonal  disks.  In  both  cases,  the 
areas  are  kept  constant  and  equal  to  that  of  a  circular  disk  with  radi¬ 
us,  A.  A  disk  monopole  radius  of  0.2A  is  used  for  the  attachment 
mode.  The  most  dominant  modes  in  this  problem  are  the  wire  and  attach¬ 
ment  modes.  This  is  especially  true  for  small  disk  radius.  As  the 
disk  radius  becomes  larger,  the  surf ace-patch  modes  become  more  influen¬ 
tial  to  the  final  result.  The  computed  results  are  compared  with  the 
data  obtained  from  the  formulation  by  Richmond  C28D.  The  result  from 
the  square  disk  approximat ion  agrees  with  Richmond's  result  up  to  about 
A  =  0.4  A  ,  while  that  of  the  octagonal  disk  approximation  is  in  good 
agreement  up  to  A  =  1.0A  .  The  formulation  used  in  this  work  can  solve 
the  problem  of  the  off-centered  monopole,  which  was  not  treated  in 
any  previous  work. 
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Figure  24.  Resistance  of  a  monopole  antenna  at  the  center  of  a  disk  in 
free  space. 


B.  Scattering  from  Polygonal  Plates 

In  this  section,  the  scattering  data  from  plates  with  different 
sizes  and  shapes  are  compared  with  independent  calculations  and  measure¬ 
ments.  The  bistatic  and  backscatter  data  are  plotted  in  terms  of  the 
radar  cross-section  (or  echo  area)  defined  as 


o=  Wm 

r  -*ao 

where  r  is  the  radial  distance  from  the  plate  to  the  observer  and  S 
and  S-  denote  the  time  average  power  densities  in  the  scattered  andsin- 
cident  fields,  respecti vely. 

The  backscatter  measurements  were  made  at  9.067  GHz.  The  plates 
were  cut  from  0.032-inch  thick  aluminum  sheet  stocks.  A  transit  was 
used  to  align  the  suspended  plate  with  the  polarization  axis  of  a  horn 
antenna.  All  these  measurements  were  done  at  The  Ohio  State  University 
ElectroScience  Laboratory  by  T.  Chu. 
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1.  Circular  Disks 

Figure  26  shows  the  broadside  radar  cross-section  of  an  octagonal 
plate  at  different  disk  radii.  Also  shown  are  the  data  obtained  from 
the  eigenfunction  solution  of  a  circular  disk  by  Hodge C29J.  Subdivi¬ 
sion  method  1  is  used  to  set  up  the  surface-patch  modes  except  at  the 
frequency  where  the  disk  radius  is  1.7  A.  Here  subdivision  2  is  used 
to  reduce  the  number  of  unknowns  and  computation  time. 

Backscatter  data  from  an  octagonal  plate  with  a  radius  of  1.2  A 
are  shown  in  Figures  27  and  28  for  the  6  and  polarization,  respect¬ 
ively.  Bistatic  data  for  the  same  plate  are  shown  in  Figures  29  and  30 
for  the  0  and  <J>  polarization,  respectively.  Comparisons  of  the  scat¬ 
tering  results  from  the  octagonal  plate  with  the  eigenfunction  solution 
by  Hodge  C29D  show  disagreements  between  the  present  calculations  and 
Hodge's  data  at  angles  close  to  and  at  grazing  incidence.  This  is 
due  to  the  fact  that  an  octagon  only  approximates  a  circle.  The  eleva¬ 
tion  plane  pattern  for  a  circular  disk  should  be  independent  of  <t>. 
However,  the  pattern  made  at  <j>  =  0°  and  <J>  =  22.5°  for  the  octagonal 
plate  are  not  identical.  Note  that  the  solution  for  the  circular  disk 
falls  between  these  two  curves.  This  seems  to  indicate  the  effect 
of  the  octagonal  approximation  to  a  circular  disk. 


2.  Triangles 

The  backscatter  data  from  a  right  angle  triangle  with  the  side 
lengths  of  3  x  are  shown  in  Figures  31  through  33.  The  surface-patch 
mode  setup  using  subdivision  method  2  and  the  Gordon-Hall  procedure  are 
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Figure  26.  Normal  incidence,  backscatter  from  disk. 


Backscatter  from  disk  with  A/>  =  1.?. 
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Figure  29.  Normal  incidence,  bistatic  scatter  from  disk  with  A/A  =  1 
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(BISTATIC)  RADAR  CROSS  SECTION,  <r/x2  ( dB) 


Figure  30.  Normal  incidence,  bistatic  scatter  from  disk  with  A/A 
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Figure  33.  Backscatter  from  a  3-a  right  angle  triangle  (e-polari¬ 
zation)  in  the  elevation  plane. 
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shown  in  Figures  34  and  35,  respectively.  Note  that  the  Gordon-Hall 
procedure  set  up  264  surface  dipole  modes,  while  subdivision  method  2 
set  up  only  176  modes.  The  backscatter  data  obtained  by  using  tie 
two-mode  layouts  and  measurements  are  shown  for  the  azimuth  plane  with 
0-  and  ^-polarization  and  the  o-pol ar ization  for  the  elevation  plane 
in  Figures  31,  32,  and  33,  respectively.  Although  the  result  from 
the  Gordon-Hall  mode  layout  is  not  shown  in  Figure  33,  it  agrees  well 
with  the  two  curves  shown. 

The  backscatter  for  an  isosceles  triangle  in  the  azimuth  plane 
for  0-  and  4>-polarization  are  shown  in  Figures  36  and  37,  respectively. 
The  surf ace-patch  mode  layout,  using  subdivision  method  2  as  shown  in 
Figure  38,  illustrates  how  the  quadri lateral  surf ace-patch  modes  are 
adequate  for  modelling  the  surface  currents  on  a  plate  with  a  very  acute 
angle.  The  acute  angle  is  30°  and  the  perpendicular  height  of  the 
triangle  is  3.25  A.  No  measurement  was  done  for  the  ^-polarization 
(see  Figure  37),  but  the  result  shown  can  be  assumed  to  be  accurate. 


3.  Fin-shaped  Plate 

A  fin-shaped  plate  along  with  the  surf ace-patch  mode  layout  is 
shown  in  Figure  39.  This  plate  geometry  again  illustrates  that  subdivi¬ 
sion  method  2  does  set  up  modes  with  relatively  uniform  sizes.  Figures 
40  and  41  show  the  calculations  and  measurements  of  the  backscatter 
in  the  azimuth  plane  for  the  0-  and  4>-polar  ization,  respectively.  In 
the  4>-polarization  (Figure  41),  a  curve  obtained  by  using  the  Geometric 
Theory  of  Diffraction  (GTD)  is  shown.  This  GTD  formulation,  which  in¬ 
corporated  equivalent  edge  diffracted  currents  C  303  ,  is  current¬ 
ly  being  developed  at  The  Ohio  State  University  ElectroScience  Labora¬ 
tory  C  313. 
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Figure  35.  A  set  of  overlapping  surf  ace-patch  inodes  formed  in  a  3-A 
right  angle  triangle  by  using  the  Gordon-Hall  procedure. 
This  set  of  modes  is  used  to  generate  the  backscatter  data 
in  Figures  31  and  3?. 


RADAR  CROSS  SECTION,  er  /  \ *  (dB) 


A  set  of  overlapping  surface-patch  dipole  modes  formed  in  a 
30°  isosceles  triangle  by  using  subdivision  method  2.  This 
set  of  modes  is  used  to  generate  the  backscatter  data  in 
Figures  36  and  37. 
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Figure  40.  Backscatter  from  a  3-X  fin-shaped  plate  (0-polarization) 
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Figure  41.  Backscatter  from  a  3-A  fin-shaped  plate  (^-polarization). 
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CHAPTER  V 


ANALYSIS  OF  MICROS  TRIP  ANTENNAS 


This  chapter  will  present  an  analysis  of  the  microstrip  antenna. 
Surface-patch  dipole  modes  are  used  to  model  the  microstrip  patch  and 
volume  polarization  currents  for  the  dielectric  slab.  The  nature  of 
this  problem  requires  unusually  precise  computation  of  the  impedance 
matrix,  but  the  method  is  capable  of  accurately  predicting  currents, 
impedance,  and  resonant  frequency  of  the  antenna. 

A.  Introduction 


A  microstrip  antenna  consists  of  a  metallic  plate  or  patch  on 
an  electrically  thin  grounded  dielectric  slab.  The  patch  shapes; 
such  as  rectangular,  circular,  semicircular,  trapezoidal,  triangular, 
and  pentagonal  r 3271 ,  have  been  fabricated  in  practical  microstrip 
antennas.  The  patch  can  be  fed  either  by  a  microstrip  transmission 
line  contacting  an  edge  of  the  patch,  or  by  a  coaxial  probe  extending 
through  the  ground  plane  and  contacting  the  patch  (see  Figures  42 
and  43).  There  can  be  more  than  one  feed  point. 

The  microstrip  antenna  lias  a  very  low  profile  since  its  height 
above  the  ground  plane  is  basically  the  dielectric  substrate  thickness. 
Except  for  the  coaxial  feedline,  the  antenna  does  not  extend  behind 
the  ground  plane.  Its  simple  construction  makes  it  very  lightweight. 
These  features  allow  microstrip  antennas  to  be  built,  conformal  to 
the  surface  of  an  air-  or  spacecraft.  Because  the  antenna  is  above  the 
ground  plane,  it  can  be  easily  mounted  to  any  vehicle  without  interfer¬ 
ing  with  its  internal  structure.  Microstrip  phased  arrays  have  also 
been  built.  In  this  case,  the  microstrip  radiators  and  their  associated 
microwave  feed  system  (microstrip  transmission  lines)  can  be  etched 
on  a  single  surface.  Thus,  these  antennas  are  inexpensive  and  simple 
to  fabricate.  Although  typical  bandwidths  are  in  the  order  of 
a  few  percent,  advantages  can  outweigh  this  shortcoming.  A  compre¬ 
hensive  report  on  the  various  aspects  ot  microstrip  antennas  is  given 
in  reference  f  32 1 . 

Since  the  early  work  on  microstrip  antennas  reported  by  Deschamps 
in  1953)  34'  up  to  1975,  the  analysis  and  design  ot  these  antennas 
wore  done  by  using  simple  models  adapted  from  transmission  line  theory. 
Howell  !  35  and  36.  modelled  the  microstrip  antenna  as  a  dielectric 
loaded  cavity  with  magnetic  side  wall.  Using  tins  model,  the  resonant 
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frequency  for  rectangular  and  circular  microstrip  patches  can  be  predict 
ed,  but  the  input  impedance  at  the  feed  has  tu  be  determined  experimental 
1 y.  8ased  on  the  same  model  and  taking  into  account  the  fringing  field, 
a  more  accurate  expression  for  the  resonant  frequency  of  a  circular  disk 
was  later  given  by  Shen,  et  a  I . ,  37  .  Munson  !  38  and  Derneryd  !  39  ! 

considered  the  microstrip  antenna  as  two  slots  perpendicular  to  the 
feedline  and  separated  by  a  very  low  impedance  transmission  line. 

The  slots  and  the  surface  patch  wore  modelled  by  an  equivalent  network, 
and  transmission  line  theory  was  used  to  calculate  the  input  impedance. 
This  method  was  used  to  model  rectangular  micros  trip  patches  fed  at 
the  center  of  an  edge.  The  theory  was  extended  by  Derneryd  f40  to 
account  for  different,  feed  locations  along  the  edge  and  the  mutual 
conductance  between  the  slots. 

More  recent  work  on  microstrip  antenna  analysis  is  based  on 
modelling  the  region  between  the  patch  and  the  ground  plane  as  a  reson¬ 
ant  cavity.  In  this  method,  a  complete  set  of  eigenvectors  which 
satisfy  the  wave  equation  and,  therefore,  can  represent  the  fields 
inside  the  cavity  must  be  found,  if  the  patch  shape  is  of  a  canonical 
form,  the  eigenvectors  can  be  written  in  terms  of  known  functions. 

In  order  to  simplify  the  boundary  condition  at  the  side  wall  (see 
Figure  43),  many  investigators  assume  that  the  tangential  component 
of  the  magnetic  field  (d)  along  the  wall  is  negligible.  Thus,  the 
side  wall  becomes  a  magnetic  wall.  Using  this  model.  Long,  et  al . , 

(411  analyzed  the  circular-disc  microstrip  antenna  and  derived  the 
expressions  for  the  far  field,  directive  gain,  efficiency,  and  the 
Q-factor  of  the  antenna.  An  extension  or  this  work  was  reported  by 
Derneryd  where  the  circular  disk  ''IP'  arid  rectangular  patch  F z) 3  were 
analyzed.  A  comprehensive  work  using  this  inode  I  was  given  by 
Lo,  et  al , ,  [441.  Here,  solutions  for  numerous  canonical  shapes  were 
given.  It  should  be  noted  that  the  fringing  field  along  the  edge 
of  the  patch  was  taken  into  account  either  by  using  an  effective  dielec¬ 
tric  constant  for  the  substrate  and/or  <-•  small  extrusion  in  the  diniens  ion 
of  the  patch.  An  alternate  approach  to  account  for  the  fringing  field 
is  the  use  of  an  impedance  boundary  condition  on  the  side  wall.  In 
other  words,  the  side  wall  is  no  longer  an  open  circuit  (magnetic 
wall),  but  has  an  associated  admittance  (or  radiating  wall).  Carver 
f  451  incorporated  this  impedance  condition  for  his  model  and  derived 
the  expression  of  the  wall  admittance  for  the  rectangular  and  circular- 
patch.  In  addition.  Carver  used  the  finite  element  method  to  solve 
the  wave  equation  for  arbitrarily  shaped  patches.  Since  the  problem 
of  finding  the  wall  admittance  for  arbitrarily  shaped  patches  was 
not  solved,  this  latter  technique  did  not  include  the  impedance  condi¬ 
tion.  C.M.  Kaloi  used  the  cavity  resonator  model  to  derive  design  equa¬ 
tions  which  can  be  used  for  a  conventional  microstrip  antenna  as  well  as 
the  "Twin  Electric  Microstrip  Antenna,"  where  identical  microstrip  patch 
es  are  formed  on  both  sides  of  a  dielectric  slab C 553. 


Agrawal  and  Bailey  ;  46  modelled  the  nnrrostrip  patch  and  feedlines 
ny  using  the  thin  wire  grid.  Using  Richmond's  thin  wire  computer 
code  [3  ,  a  solution  was  first,  obtained  for  the  antenna  radiating  in 
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the  homogeneous  medium  having  the  permittivity  of  the  dielectric  sub¬ 
strate.  The  result  was  then  shifted  in  frequency  and  the  input  imped¬ 
ance  scaled  by  constants  determined  from  measurements.  This  was  a 
moment  method  approach,  as  is  the  approach  used  in  the  following  analy¬ 
sis.  However,  surface  patch  modes  are  now  used  to  model  the  microstrip 
antenna,  and  the  presence  of  the  dielectric  is  taken  into  account 
by  the  volume  polarization  current. 

The  three  models  described  above  can  be  contrasted  by  the  types  of 
unknowns  that  have  to  be  determined.  In  the  transmission  line  model, 
the  unknowns  are  the  characteristic  impedance,  propagation  constant, 
and  loads  on  an  equivalent  transmission  line.  The  unknowns  in  the 
cavity  model  are  the  distribution  of  the  cavity  modes  (or  eigenfunc¬ 
tions)  and  the  resonant  frequencies  of  these  modes  (eigenvalues).  In 
the  moment  method  solution  presented  here,  the  unknowns  are  the  surface 
currents  flowing  on  the  microstrip  patch,  and,  if  a  dielectric  slab 
is  present,  the  equivalent  volume  polarization  currents  in  the  region 
occupied  by  the  slab. 


B.  The  Model 

In  this  section,  the  model  of  the  microstrip  patch,  the  effect 
of  the  dielectric  substrate,  and  the  numerical  difficulties  will  be 
discussed.  The  ground  plane  under  the  patch  is  assumed  to  be  infinite 
and,  thus,  the  image  theory  can  be  used.  In  this  solution,  the  equi¬ 
valent  volume  polarization  currents  in  the  dielectric  slab  are  related 
to  the  surface  charge  density  of  the  patch.  Consequently,  the  surface 
current  on  the  patch  is  the  only  unknown  involved.  As  will  be  shown 
later,  the  current  on  the  plate  is  inversely  proportional  to  the  differ¬ 
ence  of  the  self  and  mutual  impedances  of  the  patches.  Since  the 
difference  between  the  two  impedances  is  very  small,  both  impedances 
have  to  be  evaluated  very  accurately. 


1.  Air  Dielectric  Microstrip 

Figure  44  shows  the  side  view  of  an  air  dielectric  microstrip 
patch  antenna  at  a  height  T  above  a  ground  plane.  The  antenna  is 
shown  coaxially  fed,  although  it  could  just  as  easily  be  edge  fed. 
Using  image  theory,  the  ground  plane  is  removed  and  the  image  of  the 
patch  and  feed  probe  are  inserted.  The  coaxial  aperture  is  modelled 
by  a  delta  gap  voltage  generator.  This  is  shown  in  Figure  4f.  Gener¬ 
ally,  the  currents  at  the  wire/patch  junction  are  modelled  by  the 
attachment  mode  discussed  in  Chapter  Two.  Bui  our  experience,  ant' 
apparently  the  experience  of  others  Fl4  I,  indicates  tnat  the  feed 
probe  and  voltage  generator  can  be  replaced  by  an  impressed  current 
filament  given  by 

Jj  =  2  amp  (  35) 

existing  between  the  patches  and  constant  in  the  z  direction. 
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Figure  44.  Side  view  of  an  air  dielectric  microstrip  antenna  fed  by  co¬ 
axial  line. 
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Figure  45.  The  model  for  the  air  dielectric  microstrip  antenna. 


Figure  46  shows  the  geometry  that  can  be  analyzed  using  the  moment 
method  surface  patch  modelling  technique  described  in  Chapter  Two. 

The  dominant  current  on  the  patch  is  in  the  direction  where  the  dimen¬ 
sion  of  the  patch  is  about  one-half  the  wavelength  in  the  dielectric 
medium.  Thus,  only  one  current  dipole  mode  is  used  for  each  patch. 

In  the  case  where  the  patch  is  close  to  being  a  square,  the  two  ortho¬ 
gonal  polarizat ions  are  needed.  In  this  work,  this  problem  is  analyzed 
by  solving  for  each  current  polarization  independently,  and  the  total 
input  impedance  is  the  sum  of  the  individual  input  impedances.  For 
each  case,  the  impedance  matrix.  Equation  (8),  will  have  two  unknowns. 

It  is  well  known  that  the  close  proximity  of  the  plates  in  Figures 
44  through  46  will  cause  the  surface  currents  to  distribute  unevenly 
between  the  top  and  bottom  surfaces  of  the  plates.  Actually,  the 
majority  of  the  current  will  be  on  the  interior  surfaces  of  the  plates. 
This  complexity  is  avoided  in  this  present  solution  since  the  surface 
current,  in  Equation  (6),  is  the  vector  sum  of  the  currents  on 
the  top  and  bottom  surfaces  (see  Chapter  Two).  This  is  the  current 
which  radiates  and  is,  therefore,  of  primary  interest  in  any  solution. 


2.  Modification  for  Dielectric  Slab 

It  will  now  be  shown  how  the  matrix  equation  for  the  air  dielec¬ 
tric  inicrostrip  can  be  modified  to  account  for  the  dielectric  slab 
(or  substrate).  Figure  47  shows  the  geometry  of  Figure  46,  but  with 
a  dielectric  slab  of  permi tt i vi ty  c ,  or  relative  permittivity  er  - 
e/e  ,  between  the  plates.  Using  the  volume  equivalence  theorem,  the 
slat)  is  removed  and  replaced  by  free  space  and  the  equivalent  volume 
polarization  currents  (see  Figure  48) 

Jv  =  jo)(e-e0)  E  (36) 

where  E  is  the  actual  electric  field  in  the  slab.  Note  that  Jy  exists 
only  in  the  region  previously  occupied  by  the  slab.  Jy  is  unknown 
since  £  is  unknown.  Jy  can  be  treated  as  another  set  of  independent 
unknowns  and  expanded  in  terms  of  Ny  basic  functions.  The  coefficients 
’ n  the  expansion  for  and  Jy  are  determined  by  solving  a  system 
of  N  +  Ny  simultaneous  equations  C473.  However,  the  additional  un¬ 
knowns,  and  consequently  the  increase  in  computation  time,  made  this 
technique  unattractive  for  this  problem.  An  alternative  approach 
used  in  this  work  is  to  treat  Jc;  and  Jy  as  dependent  unknowns  r 48 
and  4<f).  With  the  slab,  the  RIE  (6)  becomes 

-  (J^  '  IE,||)  ds  -  J"/J  (Jy  •  Lm)  dV  ;  VMi  (  ,/) 

Vd 

where  is  the  volume  occupied  by  the  slab  and  Eni  is  still  the  field 
in  free  space.  Assuming  that  the  field  inside  the  dielectric  slab 
due  to  the  surface  patch  current  expansion  mode  can  be  found,  J  is 
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Figure  46.  The  voltage  generator  in  Figure  45  is  replaced  by  1-amp  im¬ 
pressed  current  filament. 


\ 


Figure  47.  Microstrip  antenna  with  dielectric  slab. 
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made  dependent  on  J  s by  expressing  the  currents  in  Equation  (37)  by 
the  expansion  modes 

£n=£n+£n  (38) 

where  is  the  surface  conduction  current  and  F_P  is  the  volume  polari¬ 
zation  current,  n 

In  =  >(e-c0)  Ip  (39) 

where  En  is  the  field  in  the  dielectric  due  to  the  n^1  expansion  mode, 
Fn.  Using  Equations  (36),  (38),  and  (39),  the  R1E  (37)  can  be  express¬ 
ed  by  the  matrix  equation 

N 

n  =  1  *n  (zmn  +  ^zmn^  "  %  »  m  -  2,  N  (40) 

where  7mri  and  Vm  are  defined  by  Equations  (  9)  and  (10),  respecti vely, 

.  mn  rnn 

and 


Azmn  =  -ju)(e-€0)  Iff  Ef'  •  Em  dV  (41) 

N  is  the  total  number  of  surface  patch  current  expansions.  While 
the  integral  in  Equation  (41)  in  theory  extends  over  the  entire  volume 
of  the  slab,  in  practice  it  is  sufficient  to  integrate  a  few  slab 
thicknesses  beyond  the  region  shared  by  modes  m  and  n.  Although  the 
expression  for  E^  in  terms  of  Somerfeld  integrals  can  be  found!  bO  1  , 
here  we  make  some  simplifying  approximations. 


Assuming  that  only  the  z  component  of  £p  is  significant  in  the 
slab,  the  procedure  for  computing  E^z  will  now  be  presented.  Figure 
&9  shows  _Jn  radiating  in  the  presence  of  the  dielectric  slab.  The 
surface  charge  density  associated  with  is  obtained  using  the  equation 
of  continuity 


V  -Jn  =  -ju>psr) 

(A?) 

or 

0  =  --V4n 

,Jsn  jco 

(4  3) 

where  Psn  is  the  surface  charge  density  associated  with  dr). 
the  boundary  condition 

From 

°1  +  =  8Sn 

(44) 
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49.  Jn  in  the  presence  of  the  dielectric  slab. 
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and  if  we  were  to  assume  0^  -  D-,  -  0  ,  the  flux  density  radiated  in 
a  homogeneous  medium  (regardless  of  its  permittivity),  then  Equation 
(43)  would  become 


D0  =  Psr/2  =  :1 lin  •  («5) 

2jw 

Finally,  the  electric  field  at  the  slab  surface  would  be 

Ep7  =  Do/c  '-l--'  (^) 

2  j  air 

The  major  approx imat ion  involved  in  obtaining  Equation  (46)  was  the 
assumption  0,  =  Op  =  0  .  Actually,  the  pr.'sm,  e  of  the  slab  causes 
a  significant  asymmetry;  that  is,  D,  >  D  >  Op,  Equation  (46)  can 
be  modified  to  include  this  effect  as 


Edz  =  RD  /F  =  "RV  *  — n  ( /J7) 

2  j  oae 

where  the  ratio  R  =  D,/D  can  be  evaluated  using  the  method  of  dielec¬ 
tric  images  i54  !  and  is  derived  in  Appendix  B.  Figure  50  illustrates 
a  typical  asymmetry  of  the  D  fields  at  the  surface  of  a  slab.  R 
is  evaluated  very  close  to  the  surface  patch  and  has  been  shown  to 
be  accurate  for  very  thin  slabs  (TcO.Olx  ).  However,  for  thicker 
slabs,  the  final  results,  begin  to  deteriorate  and  this  is  caused  by 
the  approximation  for  E^z. 

'"•[nation  (47)  predicts  that  Edz  at  the  surface  of  the  slab  is 
R  times  the  field  radiated  by  Jn  in  a  homogeneous  medium  of  permittivity 
r.  .  We  use  the  approx  imat  ion  that  E^z  ,  at  any  point  in  the  slab, 
is  R  times  the  field  radiated  by  Jn  in  a  homogeneous  medium  of  permit¬ 
tivity  c  . 

The  excitation  column  could  be  evaluated  using  Equations  (10) 
and  (35),  i.e. 


2T 

V m  -  J  Em  ’  j  d  k  ( n 

0 

i  f  is  taken  to  be  t  he  field  of  the  m^1  test  mode  in  the  presence 
of  tne  dielectric  slab.  We  would  like  to  use  equation  (48)  witn  _Em 
still  the  free  spai  .•  field  of  the  in  test  mode.  Using  the  results 
of  Equation  (47),  this  can  be  done  by  evaluating  the  Vm  using  Equation 
(  4g)  and  then  imi  1  lip  I  vinj  by  R/  c_ . 
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RELATIVE  D 


y 


Figure  50.  Relative  field  along  the  line  x  =  3.75  cm.  and  y 
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Once  the  patch  currents  are  known,  it  is  straight  forward  to 
evaluate  the  input  impedance  at  the  feed  point.  Since  the  feed  current 
is  taken  as  one  amp,  the  input  impedance  is  numerically  equal  to  the 
voltage  at  the  feed  point,  or 


2in 


N 

E 

n  =  1 


!n  Vn 


(49) 


where  the  factor  h  is  a  result  of  the  use  of  image  theory  to  model 
the  antenna.  If  the  antenna  is  coaxially  fed  (rather  than  edge  fed 
by  microstrip  line),  the  effects  of  the  probe  can  be  included  by  adding 
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J'xn 

'  JVTT 

tan 

to  zin. 

3. 

Expansion 

and 

The 

expression 

for 

given  by 

Equation  (12).  i 

(so) 


many,  we  choose  poeQ  =2n/AQ>  where  A 

is  the  free-space  wavelengthL  13  and  16].  In  theory,  the  moment  method  0 
should  converge  to  the  same  result  regardless  of  any  reasonable  choice 
for  k.  Here  we  choose 


k  s  “V">v7 


(51) 


where  e  is  the  effective  dielectric  constant  given  by  £51] 


.+1 


4^  (i  ♦  — i 

Z  go 


(5?) 


where  w  is  the  width  of  the  microstrip.  The  choice  of  k  primarily 
effects  the  convergence  of  the  impedance  level  rather  than  the  resonant 
frequency.  This  convergence  can  be  especially  slow  for  a  poor  choice 
of  k  since  the  impedance  level  is  essentially  proportional  to  the 
divergence  or  slope  of  the  microstrip  patch  surface  current  density 
at  the  feed  point. 


4.  Numerical  Difficulties 

The  moment  method  solution  for  microstrip  antennas  requires  the 
unusually  precise  computation  of  the  elements  in  the  impedance  matrix. 
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Consider  the  air  dielectric  inicrostrip  of  Figure  46.  If  we  let  the 
subscript  b  refer  to  a  single  mode  on  the  bottom  or  image  patch  and 
the  subscript  t  refer  to  a  single  mode  on  the  top  patch,  then 


bb  +  *t  Zbt  ^b 

(  d3  ) 

tb  +  *t  Ztt  =  Vt 

(54) 

If  we  assume  that  modes  b  and  t  are  both  pointing  to  the  right,  then 
from  image  theory  we  know  1^  =  -1^.  Also,  from  the  symmetry  of  the 
problem  Z.  .  =  Zl.  ,  Zfc.  =  Z.  t,  and  V.  =  -V.  .  In  this  case.  Equations 
(5?)  and  f 54 )  can  be  solved  as 


h  -  -Ib 


*tt  "  Ztb 


(t,E) 


Equation  (55)  shows  that  the  solution  is  critically  dependent  upon 
the  difference  between  the  self  and  mutual  impedance  of  the  top  ami 
bottom  inodes.  Since  modes  b  and  t  are  identical,  except  that  they 
are  shifted  by  typically  0.01  X,  we  expect  that  Ztt  ~  z^.  Thus, 
obtaining  Ztt  -  Ztb  accurately  is  a  difficult  numer.cal  problem. 

As  an  example,  Table  1  shows  Z^t,  Ztb  and  related  quantities 
for  the  inicrostrip  shown  in  Figure  5T,  except  that  e  =  1.  Data  are 
shown  at  the  three  frequencies  952.0,  963.0,  and  974.0  MHz  which  are 
below,  at,  and  above  resonance,  respectively.  (The  data  in  Table 
1  and  2  do  not  include  contributions  from  losses  in  the  plates  or 
the  dielectric. ) 

Note  that  the  self  resistance  is  always  slightly  larger  than 
the  mutual  resistance.  This  difference  controls  the  level  of  the 
plate  currents  and  of  the  input  impedance,  with  both  increasing  as 
the  difference  becomes  smaller.  The  self  reactance  is  smaller  tfian 
the  mutual  reactance  below  resonance,  but  larger  above.  At  resonance 
the  self  and  mutual  reactances  are  equal,  and  Z , .  -  Zt.  is  purely 
real.  Since  the  elements  in  the  voltage  vector  are  essentially  pu-'r-ly 
imaginery.  Equation  (55)  shows  that,  at  resonance,  the  currents  are 
maximum  and  essentially  purely  imaginary.  Finally,  from  Equation 
(49),  it  can  be  seen  that  at  resonance  the  impedance  is  maximum  and 
real . 


Next,  consider  the  same  antenna,  except  with  the  dielectric  slab. 
Table  2  shows  Ztt,  Zib»  AZtt,  and  AZt(  and  related  quantities  at  the 
frequencies  630.0,  6J4.0,  and  638.0  MHz  which  are  below,  at, 
and  above  resonance,  respectively.  Note  that  Re(Z.t)  is  always  larger 
than  Re(Z^. ),  but  by  only  a  few  hundredths  of  an  onm.  Also,  lm(Z.t) 
is  less  tnan  I  ( Z )  by  about  16  ohms.  Thus,  without  the  AZ  dielectric 
correction  terms  we  are  very  far  from  resonance.  In  this  example, 

AZt|.  =  -  AZ^  is  about  j.°,  ohms.  With  this  correction,  +  AZtt 
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Figure  51.  Input  impedance  of  rectangular  microstrip  antenna. 
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is  nearly  equal  to  Z  ,  +  aZ^.  near  resonance.  At  resonance  (Z^t  + 

AZtt)  -  (Z^  +  AZ t ^ c i s  purely  real  and  equal  to  about  0.03  ohms. 

Tables  1  and  2  illustrate  the  accuracy  required  to  obtain  accurate 
current  or  impedance  data.  In  Table  2,  the  self  and  mutual  resistance 
must  be  evaluated  to  an  accuracy  of  about  0.05  percent  or  better, 
in  order  that  their  difference  of  0.04  ohms  be  reasonably  accurate. 

Also,  Table  2  shows  that  a  combined  error  of  jl  ohm  in  Z,t,  Zj.  ,  and 

will  lead  to  a  shift  in  the  resonant  frequency  of  a&out  27  M Hz, 
or  about  a  four  percent  shift. 

Losses  in  the  finitely  conducting  patches  and  in  the  dielectric 
can  signif icantly  effect  the  current  and  impedance  level  and  bandwidth 
of  the  micro. trip  antenna.  The  finite  conductivity  of  the  patches, 
whose  treatment  is  analogous  to  that  of  thin  wires[19],  increases 
the  self  resistance  while  leaving  the  mutual  resistance  unchanged. 

Losses  in  the  dielectric  are  treated  by  simply  allowing  c.  to  be  complex 
in  Equation  (41).  The  result  is  to  increase  the  self  resistance  and 
decrease  the  mutual  resistance  by  equal  amounts.  Thus,  in  both  cases, 
the  losses  increase  the  difference  between  the  self  and  mutual  resistan¬ 
ces.  For  reasonably  low  loss  materials,  the  changes  will  be  a  tiny 
fraction  of  an  ohm.  However,  this  change  will  be  significant,  to 
the  current  and  input  impedance,  if  it  is  comparable  to  the  difference 
between  the  self  and  mutual  resistance  fir  the  lossless  antenna;  i.e. 
for  the  example  shown  in  Table  2,  0.04  ohms.  Mote  that  this  strong 
dependence  of  the  current  level  on  the  losses  suggests  that  it  may 
be  inaccurate  to  use  the  loss  free  currents  and  fields  in  perturbation 
technique  to  evaluate  the  loss.-,  in  the  plates  and  dielectric. 


C.  Numerical  Results 


In  this  section,  computations,  based  upon  the  results  of  the 
previous  section,  will  be  presented  and  compared  to  measurements  and 
previous  computations.  All  of  the  present  results  were  made  using 
the  same  user  oriented  computer  code.  The  computations  in  Figures 
51  through  53  were  made  with  one  unknown,  while  two  unknowns  were 
required  in  Figure  54. 

Figure  51  shows  a  coaxially  fed  rectangular  microstrip  of  length 
15  cm.  and  width  7.5  cm.  The  dielectric  slab  is  1/8  inch  thick  and 
with  relative  permittivity  =  2.56.  Figure  51  shows  good  agreement 
for  the  measured  and  computed  input  impedance  in  terms  of  impedance 
level,  resonant  frequency,  and  bandwidth  (note:  the  measurements  in 
Figures  51  and  52  were  made  by  D.H.  Schauhert  at  Harry  Diamond  Labora¬ 
tories).  Also  shown  in  Figure  51  is  the  computed  impedance  for  the 
same  antenna,  except  all  losses  are  removed  (i.e.  plates  have  peer-  i 
conductivity  and  dielectric  tan  A  -  O).  Note  the  substantial  increase 
in  impedance  level.  Figure  52  shows  a  trapezoidal  microstrip  patch, 
identical  to  that  in  Figure  61,  except  that  the  width  vnri".  n  h 
5  to  10  cm.  Note  that  the  theoretical  mode)  accurately  predicts 
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mpedance  of  trapezoidal  microstrip  antenna. 


Figure  53.  Input  impedance  of  rectangular  microstrip  antenna. 


MICROSTRIP 


mpedance  of  almost  square  microstrip  antenna. 


the  upward  shift  in  resonant  frequency.  Figure  63  compares  the  input 
impedance  of  a  rectangular  edge  fed  microstrip  of  length  7.6  cm.  and 
width  11.4  cm.  computed  hy  the  present  theory  tu  that  measured  by 
Lo,  et  al . ,  r44  ..  The  dielectric  slab  is  1/16  inch  thick  with  relative 
permittivity  e  -  2.62.  Finally,  Figure  54  shows  the  input  impedance 
for  a  nearly  square  corner  fed  microstrip  of  dimension  11.3  by  11.2 
cm.,  and  with  the  same  slab  as  Figure  63.  The  cusp  in  the  curve  is 
caused  by  Ihe  resonance  associated  with  the  11.3  cm.  side  being  at 
323  MHz,  while  that  of  the  11.2  cm.  side  is  nearly  the  same  at  32'). 5 
MHz. 


D.  Discussion 

A  surface  patch  moment  method  solut  i  on  for  the  microstrip  antenna 
has  been  presented.  The  microstrip  plates  or  patches  are  modelled 
by  equivalent  surface  currents,  while  the  dielectric  slab  is  modelled 
by  equivalent  velum-  ,»  iluri ration  currents.  Numerical  examples  illus¬ 
trate  that  the  technique  can  predict  with  good  accuracy  the  input 
impedance  of  edge  and  coaxially  fed  reel  injui-ir  and  nonrectarigu  lar 
microstrip  antennas,  I.ossr  .  in  the  pin  ind  in  the  dielectric  are 
included  and  have  a  significant,  effect  on  the  bandwidth  and  impedance 
level . 

An  interesting  feature  of  the  moment  method  solution  is  the  unusual 
precision  required  in  the  computut  infs  ■■  ■  f  el  merits  in  the  impedance 

matrix.  It  is  felt  that  this  is,  in  part,  a  consequence  of  trie  nature 
of  the  moment  method  model.  In  particular,  the  fields  and  currents 
in  a  microstrip  antenna  .me  very  nearly  those  of  a  closed  cavity. 

In  the  moment  method,  the  structure  is  modelled  by  a  number  of  plates, 
which  are  basically  open  radiators.  It  is  hoped  that  by  enforcing 
boundary  conditions  on  the  plates  that  I'ne  s  >1  o  i  m  will  "see"  the 
closed  or  cavity  nature  of  the  structure.  In  fact,  it  does;  however, 
the  description  of  the  plates  and  their  location  is  required  to  bo 
very  prec  ise.  A  second  explanat  ion  o mes  I  rum  cons i derat  ion  of  liqua¬ 
tions  (53)  and  (54).  Equation  (5°)  enforces  (in  an  approximate  sense) 
the  boundary  condition  that  the  tango’d  i  -I  elo.  tri  ;  field  vanishes 
on  the  bottom  or  image  plate.  Equation  (54)  onto*’- -s  the  boundary 
condition  that  the  tangential  electric  field  vanishes  on  the  top  plate. 
Since  the  top  and  bottom  plates  are  separated  by  typically  only  0.01k  t 
Equations  (53)  and  (64)  arc  mil  ore  ing,  to  ,1  I  irge  extent,  the  identical 
conditions.  Yet  accurately  solving  Equations  (62)  and  (54)  is  dependent 
upon  them  "seeing"  the  relatively  small  dilference  between  these  condi¬ 
tions.  The  result.,  i  1  lust.rat.od  in  ;  |ial  ,  (!•'.),  is  that  tin*  solution 

is  critically  dependent  upon  the  relatively  small  dilferences  between 
the  self  and  mutual  impedances. 

Cons  i  derail  !e  comparisons  tiave  been  .11  ad*  ■  between  the  moment,  meth.id 
-,d1  ul. ion  pr.;  ,  nt  ■  I  m>r  •  ol  t r  insm i ss  i or  1  in  m!  lidil  solutions 

for  planar  rectangular  microstrips.  Although  the  moment  method  solu¬ 
tion  is  more  complex,  it  is  not  significantly  more  accurate.  The 
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advantage  of  the  moment  method  solution  may  be  in  its  adaptability 
to  the  more  involved  problems  of  conformal  or  coupled  microstrip  anten¬ 
nas. 
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CHAPTER  VI 


CONCLUSION 


A  formulation,  based  on  the  R I E  with  electric  test  source,  for 
analyzing  electromagnetic  antenna  and  scattering  problems  has  been 
presented.  The  use  of  quadrilateral  piecewise  sinusoidal  surface  patch 
patch  modes,  thin-wire  modes,  and  the  wire/plate  attachment  modes  allows 
a  large  class  of  problems  associated  with  arbitrarily  shaped  surfaces 
to  be  analyzed.  Special  attention  has  also  been  given  to  the  problem 
of  automatically  setting  the  surface-patch  dipole  modes  for  a  general 
polygon.  The  applicability  and  versatility  of  this  formulation  has 
been  demonstrated  for  a  wide  range  of  problems,  including  monopole 
antenna  on  a  disk,  scattering  from  polygonal  plates,  and  analysis 
of  microstrip  antennas. 


A.  Topics  for  Further  Work 

Although  the  computer  code  developed  for  this  work  can  handle 
any  plate  geometry,  considerable  reduction  in  the  computation  time 
could  result  if  this  code  is  integrated  into  the  well  developed  and 
efficient  user  oriented  code  for  rectangular  patches.  This  savings 
in  time  is  obtained  by  the  new  code's  ability  to  identify  the  rectangu¬ 
lar  plates  and  applying  the  topelitz-like  properties  of  the  rectangular 
surface-patches  in  evaluating  the  mutual  impedances. 

Another  topic  for  further  investigation  is  the  use  of  Sommerfeld’s 
method  to  find  the  Green's  function  of  a  surface-patch  on  a  dielectric 
substrate.  This  formulation  is  more  general  than  the  method  presented 
in  the  analysis  of  the  microstrip  antenna  in  that  it  can  treat  electri¬ 
cally  thick,  as  well  as  thin,  dielectric  slabs  and  the  patch 
can  be  on  the  slab  or  imbedded  in  the  slab. 
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APPENDIX  A 


EFFICIENT  EVALUATION  OF  PWS  FI LAMENT- TO-F I  LAMENT  IMPEDANCE 


The  mutual  impedance  between  filamentary  monopoles  with  PWS  cur¬ 
rents  can  be  evaluated  from  the  rigorous  closed  form  expression  L2 21 
or  numerical  integration.  Reference Q 22 3  refers  to  subroutine  GGMM 
which  can  be  used  to  compute  the  mutual  impedance  of  general  skew 
PWS  filamentary  monopoles  with  complex  propagation  constant, y  ,  (i.e. 
lossy  homogenous  medium).  If  free  space  is  assumed,  then  y  becomes 
pure  imaginary  and  some  of  the  exponential  integrals  in  expression 
for  the  mutual  impedance  reduce  to  cosine  integrals  which  are  simpler 
to  evaluate.  A  subroutine  with  this  simplification C52U  is  used  in 
this  work  when  the  exact  mutual  impedance  is  needed.  Since  the  near 
field  expression  from  a  PWS  filament  is  known,  only  a  one-dimensional 
integration  (on  the  expansion  filament)  is  required  to  evaluate  the 
mutual  impedance.  For  parallel  filaments,  the  computation  time  requir¬ 
ed  for  the  exact  evaluation  and  numerical  integration  (INT=2)  is  about 
the  same,  but  the  letter  is  not  accurate  at  very  small  filament  separa¬ 
tion.  Although  numerical  integration  is  also  less  accurate  than  the 
exact  evaluation  when  two  nonparallel  filaments  are  close  to  each  other, 
it  is  much  faster.  For  the  nonparallel  case,  a  five-point  ( I  NT =4 ) 
and  three-point  (INT=2)  Simpson  rule  integration  is  about  twice  and 
three  times  faster  than  the  exact  evaluation  (INT=0),  respectively. 

A  simple  test  to  decide  on  the  type  of  integration  (for  the  non¬ 
parallel  case)  is  to  find  the  distance  between  the  midpoints  of  the 
two  filaments  (Dm)  and  use  the  following  decision  role: 

Dm  >  0.25X  use  INT=2 

Dm  <  0.25 A  use  INT=0. 

However,  even  within  the  Dm  <  0.25  A  radius  there  are  many  filament 
orientations  which  require  only  INT=2.  The  following  describes  a  more 
stringent  test  which  can  substantially  reduce  the  overall  computation 
time.  Figure  55  shows  the  reactances  and  computation  times  using  INT=0, 
INT=2  and  different  INT's  chosen  by  the  following  decision  rule. 

Figure  56  shows  two  filamentary  monopoles,  M  and  N,  with  endpoints 
1,  2  and  A,  B,  respectively.  Assuming  that  M  and  N  are  not  parallel, 
the  minimum  distance,  D  -  ,  between  the  two  filaments  are  determined  by 
one  of  the  following: 


REACTANCE  (OHMS) 


Figure  55.  Reactance  of  two  PWS  fiJamentary  monopoles.  The  lengths  of 
the  filaments  are  0.25X.  SI  and  S2  denote  the  terminal  and 
end  points  of  the  test  monopole,  respectively.  T1  and  T2 
denote  the  terminal  and  end  points  of  the  expansion  monopole 
respecti vely. 


Figure  56.  The  geometry  used  to  find  the  smallest  distance  between  two 
filaments. 
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1 


A  A  A  A 

a.  If  A. «  <  0  and  ft,,  *  ft,?  <  0,  then  the  filaments  cross 

each  other.  D  .  is  trie  perpendicular  distance  between  the 

two  parallel  planes  formed  by  the  two  filaments. 

A  A  A  A 

b.  If  pMA  *  PMg  <  0  and  pN1  *  pn2  >  0,  then  Dm1n  is  the  minimum 

value  of  (|pn1| ,  |pN2I ). 


c. 


d. 


If  PNi  ‘  P^O  and  £MA  *  PMB  1  0,  then  Dmin  is  the  minimum 
value  of  (|pmaI,  I PMB I ) - 


If  PMA  "  PMB  -  0  andP 
of  tne  distances  from 


2  to  8. 


>  0,  then  D  :  is  the 
to  A,  1  to  8,2  to  A, 


minimum 

and 


D  .  is  then  used  in  the  following  decision  rule: 


0  <  D  .  <  0.03A 

—  mm  — 

choose 

INT=0 

0.03A< 

Dmini°-^ 

choose 

INT=4 

D  •  > 

min 

0.1X 

choose 

INT=2 
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APPENDIX  B 


EVALUATION  OF  R  BY  THE  METHOD  OF  DIELECTRIC  IMAGES 


The  assymetry  of  D,  and  Dp  (see  Figures  49  and  50)  due  to  the 
dielectric  slab  is  accounted  for  by  a  constant  R.  This  constant  is 
the  ratio  of  two  average  values: 

ru  C  ■ E™  <■* 

R  _  _  (56) 

;;;  ez  •  ez  dv 

J  n  m 

V2 

where  E^z  and  Ez  are  the  fields  of  the  nth  expansion  mode  evaluated  in 
the  presence  ofnthe  dielectric  slab  and  free-space,  respectively.  Ez 
is  the  field  from  the  rrrn  test  mode  evaluated  in  free  space.  V~  is  the 
common  volume  shared  by  the  modes  m  and  n.  For  electrically  thin  di¬ 
electric  slab,  the  fields  are  constant  along  the  z-coordinate  and  the 
integrals  in  Equation  (56)  reduce  to  surface  integrals.  In  practice, 
the  surface  integration  is  done  inside  Vp  and  close  to  the  surface 
of  the  nrn  mode.  Note  that  the  integrals  in  Equation  (56)  are  the  same 
as  those  used  to  determine  Azmn  in  Equation  (41). 


Since  R  is  a  ratio  of  the  field  averages,  the  static  approximation 
is  used  to  calculate  these  fields.  Consider  a  finite  line  charge, 
as  shown  in  Figure  57,  with  a  uniform  charge  distribution  of  density 
Q/h.  The  potential  distribution  of  this  line  charge  i s  C 53 ,  p.  1073 


tm4r  *  [ln  (^). 

where  =  V (x-x1)2+(y-y1)2+(z+2nd)2 
r2n  =  V (x-x2)2+(y-y2)2+(z+2nd)2 


(57) 


n  are  the  intergers  (...  -2,  -1,  0,  1,  2,  ...) 
d  is  the  dielectric  thickness  (2T) 

U  is  the  distance  from  (xp,  y?  z2)  to  the  projection  of  (x, 
y,  z)  on  to  the  line  formed  by  extending  the  line  charge,  i . e. : 
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U  =  (x-x2)  *  CAS  +  (y-y2)  *  CBS  +  (z-z2)  *  CGS 

CAS  =  (xj-x2)/h 
CBS  =  (yry2)/h 
CGS  =  (zj-z2)/h. 

Ez  is  found  by: 


% 

-1  2n 

Wg  h 


z+2nd 

.^n 


z+2nd  _ 


Equation  (58)  is  the  expression  for  Ez  from  a  line  charge  and  can  be 
integrated  to  find  the  field  of  a  surface  charge  Ez  and  Ez.  Note  that 
Equation  (56)  is  used  only  when  V2  is  the  common  area  shared  only  by 
a  monopole  from  each  mode  m  and  nf  If  both  monopoles  of  modes  m  and  n 
overlap,  then  Ez  must  be  the  field  from  a  dipole  with  +Q  on  one  monopole 
and  -Q  on  the  other.  The  field  from  a  line  charge  -Q  can  be  found 
by  the  same  method  outlined  above  and  this  is  added  to  Equation  (58) 
to  form  the  field  from  a  dipole.  The  fields  of  the  line  charges  are 
integrated  to  find  the  field  due  to  a  surface  charge. 


E2  due  to  a  surface  charge  radiating  in  a  homogenous  medium  with 
permittivity  e  has  bean  derived.  The  method  of  dielectric  images  C54J 
is  now  used  to  find  E°  ,  the  field  in  the  presence  of  a  dielectric  slab. 
Figure  58  shows  the  images  of  a  charge  Q  close  to  a  dielectric  slab 
and  the  regions  in  which  the  equivalent  images  are  applicable.  That 
is,  if  the  field  point  is  in  a  region,  then  the  fields  from  the  images 
(applicable  to  that  particular  region)  are  added.  This  technique  was 
also  used  to  generate  the  relative  Dz  in  Figure  50. 
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DIELECTRIC  SLAB, 
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Figure  58.  The  images  of  a  charge,  Q,  located  near  a  dielectric  slab. 


APPENDIX  C 


BASIC  APPROACHES  USED  IN  SUBDIVISION  METHODS  1  AND  2 


The  rectangle  is  the  simplest  shape  that  can  be  subdivided  into 
rectangular  surface-patch  modes  which  satisfy  all  the  criteria  given 
in  Chapter  Three.  Determination  of  the  set  of  modes  for  a  rectangular 
plate  like  the  one  shown  in  Figure  59  is  the  motivation  for  the  gener¬ 
al  approach  used  in  subdivision  methods  1  and  2.  The  rectangular 
is  first  divided  into  long  strips  (Figure  60)  and  these  strips  are 
then  sectioned  to  form  rectangular  blocks  (Figure  61).  These  small 
blocks  become  the  monopoles  of  the  overlapping  PWS  surf ace-patch 
modes. 

A  general  polygon  can  be  subdivided  in  a  similar  manner.  The 
polygon  is  first  divided  into  trapezoidal  strips.  These  strips  are 
then  divided  into  blocks.  It  is  the  strategies  used  to  subdivide  the 
strips  in  which  method  I  differs  from  method  2.  In  method  1,  all  the 
strips  are  divided  into  the  same  number  of  blocks  and,  thus,  generate 
a  grid  for  the  two  current  polarizations  as  shown  in  Figures  16  and 
18.  Realizing  that  the  shorter  strips  do  not  require  as  many  blocks 
as  the  longer  strips,  method  2  segments  each  strip  by  a  minimal  number 
of  blocks.  Also,  the  sides  of  the  blocks  that  do  no  coincide  with 
the  plate  boundary  are  manipulated  so  that  these  blocks  can  become 
rectangular  blocks.  In  this  way,  subdivision  method  2  can  generate  a 
large  percent  of  rectangular  modes  and  this  becomes  more  evident  in 
electrically  large  plates.  Examples  of  mode  layout  using  subdivision 
method  2  are  shown  in  Figures  23,  34,  38,  and  39.  Both  subdivision 
schemes  are  discussed  in  more  detail  in  Chapter  Three. 
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